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ABSTRACT 

CO line emission represents the most accessible and widely used tracer of the molecular interstellar 
medium. This renders the translation of observed CO intensity into total H2 gas mass critical to 
understand star formation and the interstellar medium in our Galaxy and beyond. We review the 
theoretical underpinning, techniques, and results of efforts to estimate this CO-to-H2 "conversion 
factor," A"coi in different environments. In the Milky Way disk, we recommend a conversion factor 
Xqo = 2 X 10^*^ cm~^ (K km s"'^)""'^ with ±30% uncertainty. Studies of other "normal galaxies" return 
similar values in Milky Way-like disks, but with greater scatter and systematic uncertainty. Departures 
from this Galactic conversion factor are both observed and expected. Dust-based determinations, 
theoretical arguments, and scaling relations all suggest that ATco increases with decreasing metallicity, 
turning up sharply below metallicity « 1/3-1/2 solar in a manner consistent with model predictions 
that identify shielding as a key parameter. Based on spectral line modeling and dust observations, 
Xqo appears to drop in the central, bright regions of some but not all galaxies, often coincident with 
regions of bright CO emission and high stellar surface density. This lower Xqq is also present in the 
overwhelmingly molecular interstellar medium of starburst galaxies, where several lines of evidence 
point to a lower CO-to-Il2 conversion factor. At high redshift, direct evidence regarding the conversion 
factor remains scarce; we review what is known based on dynamical modeling and other arguments. 
Subject headings: ISM: general — ISM: molecules — galaxies: ISM — radio lines: ISM 
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1. INTRODUCTION 

Molecular hydrogen, II2, is the most abundant 
molecule in the universe. With the possible exception 
of the very first generations of stars, star formation is fu- 
eled by molecular gas. Consequently, H2 plays a central 
role in the evolution of galaxies and stellar sys tems (see 
the recent review bv lKennicutt fc Evans! I201I) . Unfor- 
tunately for astronomers interested in the study of the 
molecular interstellar medium (ISM), cold II2 is not di- 
rectly observable in emission. H2 is a diatomic molecule 
with identical nucleii and therefore possesses no perma- 
nent dipole moment and no corresponding dipolar rota- 
tional transitions. The lowest energy transitions of H2 
are its purely rotational quadrupole transitions in the 
far infrared at A = 28.22 /im and shorter wavelengths. 
These are weak owing to their long spontaneous decay 
lifetimes Tdecay ^ 100 years. More importantly, the two 
lowest para and ortho transitions have upper level ener- 
gies E /k ~ 510 K and 1015 K above ground ()Dabrowskil 
|1984[) . They are thus only excited in gas with T > 100 K. 
The lowest vibrational transition of H2 is even more dif- 
ficult to excite, with a wavelength A — 2.22 /im and 
a corresponding energy E/k = 6471 K. Thus the cold 
molecular hydrogen that makes up most of the molecu- 
lar ISM in galaxies is, for all practical purposes, invisible 
in emission. 

Fortunately, molecular gas is not pure H2. Helium, be- 
ing monoatomic, suffers from similar observability prob- 
lems in cold clouds, but the molecular ISM also contains 
heavier elements at the level of a few x 10"'' per H nu- 
cleon. The most abundant of these are oxygen and car- 



bon, which combine to form CO under the conditions 
prevalent in molecular clouds. CO has a weak perma- 
nent dipole moment (/x « 0.11 D = 0.11 x 10~^^ esu cm) 
and a ground rotational transition with a low excitation 
energy hv/k « 5.53 K. With this low energy and crit- 
ical density (further reduced by radiative trapping due 
to its high optical depth), CO is easily excited even in 
cold molecular clouds. At a wavelength of 2.6 mm, the 
J = 1 — transition of CO falls in a fairly transparent 
atmospheric window. It has thus become the workhorse 
tracer of the bulk distribution of H2 in our Galaxy and 
beyond. 

As a consequence, astronomers frequently employ CO 
emission to measure molecular gas masses. The standard 
methodology posits a simple relationship between the ob- 
served CO intensity and the column density of molecular 
gas, such that 



N(H2) = Xco W(^^C^''0 J 1 ^ 0), 



(1) 



where the column density, N(H2), is in cm~^ and the 
integrated line intensity, W(CO)EI, is in traditional radio 
astronomy observational units of Kkms^^. A corollary 
of this relation arises from integrating over the emitting 
area and correcting by the mass contribution of heavier 
elements mixed in with the molecular gas. 
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^ Heneeforth we refer to the most common ^^C^^'O isotopologue 
as simply CO, and unless otherwise noted to the ground rotational 
transition J = 1 — > 0. 
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Here Mmoi has units of Mq and Lqq is usually ex- 
pressed in K km s~^ pc^. Lqo relates to the ob- 
served integrated flux density in galaxies via Leo = 
2453 S'coAu Dl/{1 + z), where S'co Au is the integrated 
line flux density, in Jy km s^^, Z?l is the luminosity dis- 
tance to the source in Mpc, an d z is the redshift (e.g., 
ISolomon fc Vanden Bout]l2005l use Eq. [7] to convert 
between W(CO) and S'coAw). Thus aco is simply a 
mass-to-light ratio. The correction for the contribution 
of heavy elements by mass reflects chiefly helium and 
amounts to a ~ 36% correction based on cosmological 
abundances. 

Both Xco and aco are referred to as the "CO- 
to-H2 conversion factor." For Xco = 2 x 
10^° cm~^(K km s~^)~^ the corresponding aco is 
4.3 M0 (K km s~^ pc^)~^. To translate integrated flux 
density directly to molecular mass, Equation [2] can be 
written as 



1.05 X 10* 



CO 



2 X 1020 



K km s 



ScoAvDj 

(1 + ^) 



For convenience we define 



Xco 



1 X 1020cm-2(K km s" 



(3) 



(4) 



We discuss the theoretical underpinnings of these equa- 
tions in S21 

Note that the emission from COJ=l^Gis found 
to be consistently optically thick except along very low 
column density lines-of-sight, as indicated by ratios of 
^^CO to ^'^CO intensities much lower than the isotopic 
ratio. The reason for this is simple to illustrate. The 
optical depth of a CO rotational transition is 



3n gj \ / Av 



(5) 



where J and Nj are the rotational quantum number and 
the column density in the upper level of the J — >■ J — 1 
transition, v is the frequency, T^x is the excitation tem- 
perature (in general a function of J, and restricted to be 
between the gas kinetic temperature and that of the Cos- 
mic Microwave Background), A w is the velocity width, /i 
is the dipole moment, gj = 2J+1 is the statistical weight 
of level J, and h and k are the Planck and Boltzmann 
constants respectively. Under typical conditions at the 
molecular boundary, r « 1 for the J ^ 1 — > transition 
requires N(H2) w 2 — 3 x lO^" cm"^ for a Galactic carbon 
gas-phase abundance Ac ~ 1.6x 10"'* (Sofia et al.:2004). 
At the outer edge of a cloud the carbon is mainly C+, 
which then recombines with electrons to form neutral C 
(Fig. [T]). Carbon is converted to CO by a series of re- 
actio ns initiated by the cosmic-ra y ionization of H or H2 
(e.g. . Ivan Dishoeck fc Blac3ll988[ ) and becomes the dom- 
inant carrier of carbon at Ay ~ 1 — 2. The CO J = 1 — )■ 
line turns optically thick very quickly after CO becomes 
a significant carbon reservoir, over a region of thickness 
AAv ^ 0.2 — 0.3 for a typical Galactic dust-to-gas ratio 

(c.f., Eq. ini). 

Equations[T]and[2]represent highly idealized, simplified 
relations where all the effects of environment, geometry. 




Fig. 1. — Calculated cloud structure as a function of optical depth 
into the cloud. Top panel shows the fractional abundance of HI, 
H2, C"*", C, and CO. Middle panel shows their integrated column 
densities from the cloud edge. Bottom panel shows the emergent 
line intensity in units of K km for [CII] 158 /xm, [CI] 609 
fim, and CO J = 1 — > 0. The grey vertical bar shows where CO 
J = 1 — > becomes optically thick. At the outer edge of the cloud 
gas is mainly Hi. H2 forms at Ay ~ 0.5 while the carbon is mainly 
C+. The C"*" is converted to C at Ay ~ 1 and CO dominates at 
Ay > 2. The model uses constant H density n = 3 X 10^ cm~^, 
a radiation fi eld x = 30 times the interstellar radiation field of 
IDrainel | |1978I '|. a primary cosmic-ray ionization rate of 2 X 10~^^ 
s~^ per hydrogen nucleon, and are base d on the PDR models o f 
IWolfire. HoUenbach fc McKeel 1 I20T0I) and lHoUenbach et al.l ||20T^ '). 



excitation, and dynamics are subsumed into the Xco 
or aco coefficients. A particular example is the effect 
that spatial scales have on the CO-to-H2 conversion fac- 
tor. Indeed, for the reasons discussed in the previous 
paragraph, Xco along a line-of-sight through a dense 
molecular cloud where Ay > 10 is not expected to be the 
same as Xco along a diffuse line-of-sight sa mpling mostly 
mate r ial where Ay < 1 (see, for example. IPineda et al.l 
120101 : iLiszt fc Petvll2012( ). Thus on smah spatial scales 
we expect to see a large variability in the CO-to-H2 con- 
version factor. This variability will average out on the 
large spatial scales, to a typical value corresponding to 
the dominant environment. Because of the large optical 
depth of the CO J = 1 — )■ transition the velocity disper- 
sion giving rise to the width of the CO line will also play 
an important role on W(CO), and indirectly on Xco- 
Indeed, there is not one value of Xco that is correct and 
applicable to each and every situation, although there are 
values with reasonable uncertainties that are applicable 
over large galactic scales. 

The plan of this work is as follows: in the remainder of 
this section we provide a brief historical introduction. In 
^we present the theoretical background to the CO-to- 
H2 conversion factor. In 21 review the methodology 
and measurements of Xco in the Milky Way, the best 
understood environment. Wc characterize the range of 
values found and the underlying physics for each mea- 
surement technique. In ^ we review the literature on 
Xco determinations in "normal" star-forming galaxies 
and discuss the techniques available to estimate Xco 
in extragalactic systems. In ^ we consider the effect 
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of metallicity, a key local physical parameter. In JJ7] we 
review the nieasurements and the physical mechanisms 
affecting the value of the CO-to-H2 conversion factor in 
the starburst environments of luminous and ultralumi- 
nous galaxies. In fJ5]we consider the explicit case of ^co 
in high redshift systems, where a much more restricted 
range of observations exist. In ^J3]we discuss the results 
of recent calculations of molecular clouds including the 
effects of turbulence and chemistry. Finally, in ^we will 
offer some recommendations and caveats as to the best 
values of ^co to use in different environments, as well 
as some suggestions about open avenues of research on 
the topic. 

1.1. Brief Historical Perspective 

Carbon monoxide was one of the first interstel- 
lar mediu m molecules observed at millime ter wave- 
lengths. iWilson. Jefferts fc PenziasI ()1970t ) reported 
the discovery of intense CO emission from the Orion 
nebula using the 36 foot NRAO antenna at Kitt 
Peak, Arizon a. Surveys of molecular clouds in the 
Galaxy (e .g., Solomon ct al. 1972; 'Wilson et aH 119741 : 
iScoville "fc Solomon 1975; Burton ct al. 1975) established 
molecular gas to be widespread in the inner Milky Way 
with a distribution that resembles giant HII regions more 
closely than that of atomic hydrogen gas. The com- 
bination of CO and 7-ray observations demonstrated 
that H2 dominates over HI by mass in the inner Galaxy 
(jStecker et al.l ll975'). By the end of the following decade, 
these studies e xtended to comple te the mapping of the 
Galactic plane ()Dame et al.lll987t) . 

The first extragalactic detections of C O occurred in 
parallel with these early Galacti c surveys (jRickard et al.l 
[1975I: ISolomon fc de ZafrallT975h . They found CO to be 
particularly bright in galaxies with nuclear activity such 
as M 82 and NGC 253. The number of extragalactic 
CO observations grew rapidly to inclu de several hundred 
galaxies over the next two decades ([Young fc Scovilld 
Il99lt I Young "eFaI1[T995h . and CO emission was employed 
to de termine galaxy molecular masses ([Young fc Scoville 
Il982t ). By the late 1980s, the first millimeter interferom- 
eter s spatially resolved molecular c l ouds in other galax- 
ies (IVogel. Boulanger fc Balll 119871: IWilson et al.lll98l . 
Such observations remain challenging, though powerful 
new interferometric facilities such as the Atacama Large 
Millimeter Array (ALMA) will change that. 

The first detection of CO at cosmological red- 
shifts targeted ultraluminous infrared sources 
and revealed very lar ge reservoirs of high l y ex- 
cited molecular gas (iBrown fc Vanden BoutI 119911 : 
ISolomon. Downes fc Radfordlll992D . Because of the deep 
integrations required, the nu mber of high redshift CO 
detections grew slowly at first ( Solomon fc Vanden Bou^ 
[2n05h . but this field is now developing rapidly driven 
by recent improvements in telescope sensitivity (see the 
review by Carilli fc Walter in this issue). An increased 
appreciation of the roles of gas and star formation in 
the field of galaxy evolution and the according need 
to determine accurate gas masses provides one of the 
motivations for this review. 

1.2. CO Excitation 

Under average molecular cloud conditions, CO 
molecules are excited through a combination of colli- 



sions with II2 and radiative trapping. They de-excite 
through spontaneous emission and collisions, except at 
very high densities where collisions are extremely fre- 
quent. Neglecting the effect of radiative trapping, ra- 
diative and coUisional de-excitation will balance for a 
critical density ricrj = Aj /^j{Tkin) (neglecting the ef- 
fects of stimulated emission), where Tkin is the kinetic 
gas temperature. Thus, for n ^ f^cr,J and excitation 
temperatures Tex,j > Ej/k « 5.53 J( J -I- l)/2 K the 
upper level of the J — ^ J — 1 transition will be popu- 
lated and the molecule will emit brightly. In these ex- 
pressions Aj is the Einstein coefficient for spontaneous 
emission (only transitions with |AJ| = 1 are allowed), 
Aj = 647rV5/iV/(3/ic35j) {Ai « 7.11 x 10"** s-^). The 
parameter "/j(T) is the corresponding coUisional coeffi- 
cient (the sum of all coUisional rate coefficients for tran- 
sitions with upper level J), which is a weak function of 
temperature. For CO, 71 ^ 3.2 6 x 10~^^ cni^ s~^ for col- 
hsions with H2 at Tk^n « 30 K (|Yang et al.ll2010[ ). Te^j 
refers to the excitation temperature, defined as the tem- 
perature needed to recover the relative populations of the 
J and J — 1 levels from the Boltzmann distribution. In 
general, T^xj will be different for different transitions. 

The critical density for the CO J = 1 transition 
is n(Il2)cr,i ~ 2200 cm~^. Higher transitions require 
rapidly increasing densities and temperatures to be ex- 
cited, as ricr.j oc and Ej cx J^. The high optical depth 
of the CO emission relaxes these density requirements, as 
radiative trapping reduces the effective density required 
for excitation by a factor ~ 1 /tj (the precise factor cor- 
responds to an escape probability and is dependent on 
geometry) . 

The Rayleigh- Jeans brightness temperature, Tj, mea- 
sured by a radio telescope for the J — > J — 1 transition 
will be 

Tj^5.53J{l~e-^') ) K. 

\e^=-.' - 1 e2-^3(.+i) _ ly 

(6) 

The final term accounts for the effect of the Cosmic Mi- 
crowave Background at the redshift, z, of interest. Note 
that frequently the Rayleigh- Jeans brightness tempera- 
ture is referred to as the radiation temperature. Ob- 
servations with single-dish telescopes usually yield an- 
tenna temperatures corrected by atmospheric attenu- 
ation (rjj), or main beam temperatures (Tmb), such 
that Tmb = Vmb where tjmb is the main beam effi- 
cien cy of the telescope at th e frequency of the observation 
(see IKutner fc UlichI I1981L for further discussion) . The 
Rayleigh- Jeans brightness temperature Tj is identical to 
Tmb for compact sources, while extended sources may 
couple to the antenna with a slightly different efficiency. 

Extragalactic results, and measurements with interfer- 
ometers, are frequently reported as flux densities rather 
than brightness temperatures. The relations between 
flux density (in Jy) and Rayleigh- Jeans brightness tem- 
perature (in K) in general, and for CO lines, are 

Sj ^ 73.5 X 10"^ 6*2 Tj w 10.9 x 10^^ 9'^ T, , (7) 

where 9 is the half-maximum at full width of the tele- 
scope beam (in arcsec), and A is the wavelength of a 
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transition (in mm). 

The Rayleigh- Jeans brightness temperature, Tj, ex- 
citation temperature, Tex,j, and kinetic temperature, 
Tkin, are distinct but related. The ^^CO transitions usu- 
ally have Tj ^ 1, making the Rayleigh- Jeans bright- 
ness temperature a probe of the excitation temperature, 
Tj ~ Tex, J for Tex J ^ 5.53 J K. In general T^x^j can be 
shown to be in the range Tcmt ^ Tex,j < T/cm- At densi- 
ties much higher than ncr,j the population of the levels J 
and lower will approach a Boltzmann distribution, and 
become "thermalized" at the gas kinetic temperature, 
Tex, J ~ Tkin- The corresponding Rayleigh- Jeans bright- 
ness can be computed using Eq. [Bl When Tex,j < Tkin, 
usually Tex, J /Tex. J ~i < 1 for lines arising in the same 
parcel of gas, and the excitation of the J level is "subther- 
mal" (note that this is not equivalent to Tj/Tj-i < 1, 
as is sometimes used in the literature). 

2. THEORETICAL BASIS 

At its core the Xco factor represents a valiant effort to 
use the bright but optically thick transition of a molecu- 
lar gas impurity to measure total molecular gas masses. 
How and why does Xqq work? 

2.1. Giant Molecular Clouds 

Because the ^^CO J = 1 — ?> transition is generally op- 
tically thick, its brightness temperature is related to the 
temperature of the rco — 1 surface, not the column den- 
sity of the gas. Information about the mass of a self- 
gravitating entity, such as a molecular cloud, is conveyed 
by its line width, which reflects the velocity dispersion of 
the emitting gas. 

A simple and exact argument can be laid for virialized 
molecular clouds, that is clouds where twice the inter- 
nal kinetic energy equ als the potential energy. Following 
iSolomon et al.l (jl987f ). the virial mass Myir of a giant 
molecular cloud (GMC) in M© is 



3(5 -2fc) , 



(8) 



where R is the projected radius (in pc), a is the ID 
velocity dispersion (in km s""'^; ctsd = V^a), G is the 
gravitational constant (G ~ 1/232 Mg^ pc km^ s~^), 
and k is the power-law index of the spherical volume 
density distribution, p{r) oc r~'^ . The coefficient in 
front of Ra^ is only weakly dependent on the density 
profile of the virialized cloud, and corresponds to ap- 
proximately 1160, 1040, and 700 for A; = 0, 1, and 2 
respectivelv (iMacLaren Richardson fc Wolfendale|[T988l : 
iBertoldi fc McKeelll992[ ). Unless otherwise specified, we 
adopt k = I for the remainder of the discussion. This 
expression of the virial mass is fairly robust if other 
terms in the virial theo rem (jMcKee fc Zweibell 119921 : 
lBallesteros-Paredesl2006[ ) . such as magnetic support, can 
be neglected (for a more general expression applicable 
to s pheroidal clouds and a g eneral density distribution 
see IBertoldi fc McKeill992l ). As long as molecular gas 
is dominating the mass enclosed in the cloud radius and 
the cloud is approximately virialized, M^ir will be a good 
measure of the H2 mass. 

Empirically, molecula r clouds are observed to follow a 
size- line width relation ()Larsonlfl981l : iHever et al.l 120091) 
such that approximately 



(9) 



with C ^ 7 km s~^ pc~° '^ (ISolomon et all 119871 : 
ISc oville et all 119871: IRoman-Duval et al.ll2010D . This re- 
lation is an expression of the equilibrium supersonic tur- 
bulence conditions in a highly compressible medium, 
and it is though t to apply under very g eneral condi- 
tions (see §2.1 of iMcKee fc Ostrikeil 120071 . for a discus- 
sion). In fact, to within our current ability to measure 
these two quantities such a relation is also approximately 
followed by ext ragalactic GMCs in galaxy disks (c 
Rubio. Leoueux fc Boulanger 1993: .Bolatto et al.. .200 



Hughes et al.ll2010D . 



Note that insofar as the size dependence of Eq. [5] is 
close to a square root, the combination of Eqs. |8]and[9] 
yields that Myir oc a*, and molecular clouds that fulfill 
both relations have a characteristic mean surface density, 
SgmCj at a value related to the coefficient of Eq. [9] so 
that Sqmc = Myir/iiR^ ~ 33 IC^ for our chosen density 
profile p cx r~^. We return to the question of Sgmc in 
the Milky Way in g31 

Since the CO luminosity of a cloud, Lqq, is the prod- 
uct of its area (tt R^ ) and its integrated surface brightness 

[Tb V2^cr), then Leo = V^^Tb <7 R^ , where Tb is the 
Rayleigh- Jeans brightness temperature of the emission 
(see Vl.21 . Using the size-line width relation (Equation 
O to substitute for R implies that Leo oc Tbc^ ■ Em- 



ploying this relation to replace a in 
a relation between M^r and Lcoj 



cx a we obtain 



mol 



200 



CO 



Tb 



0.8 



(10) 



That is, for GMCs near virial equilibrium with approxi- 
mately constant brightness temperature, Tb, we expect 
an almost linear relation between virial mass and lumi- 
nosity. The numerical coefficient in Equation [TU] is only 
a weak function of the density profile of the cloud. Then 
using the relation between C and Sqmc we obtain the 
following expression for the conversion factor 



aco 



6.1Lco"°'T, 



^co 



-0.8 ■ 



^GMC 



0.6 



(11) 



Equal ions 110! and II 1! rest on a number of assumptions. 
We assume 1) virialized clouds with, 2) masses domi- 
nated by H2 that 3) follow the size-line width relation 
and 4) have approximately constant temperature. Equa- 
tion [TT] applies to a single, spatially resolved cloud, as 
Sgmc is the resolved surface density. 

We defer the discussion of the applicability of the virial 
theorem to §4.1.1} and the effect of other mass compo- 
nents to §2.31 The assumption of a size-line width rela- 
tion relies on our understanding of the properties of tur- 
bulence in the interstellar medium. The result a oc 
follows our expectations for a highly compressible tur- 
bulent fiow, with a turbulence injection scale at least 
comparable to GMC sizes. The existence of a narrow 
range of proportionality coefficients, corresponding to a 
small interval of GMC average surface densities, is less 
wel l un derstood (for an alternative view on this point, 
see iBallestcros-Paredes et al..,2011i) . In fact, this narrow 
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range could be an artifact of th e small dynamic range of 
the samples (jHever et al.ll2009D. Based on observations 
of the Galactic center jOka et al."2001') and starburst 
galaxies (e.g., iRosolows kv & Blitz 2005), Sgmc likely 
does vary with environment. Equation [TT] implies that 
any such systematic changes in Sqmc will also lead to 
systematic changes in XqOj though in actual starburst 
environments the picture is more complex than implied 
by Equation [TTJ Section [7] reviews the case of bright, 
dense starbursts in detail. 

This calculation also implies a dependence of Xco on 
the physical conditions in the GMC, density and tem- 
perature. Combining Eqs. [S] and [S] with M^r oc pR^, 
where p is the gas density, yields a oc Meanwhile, 
because CO emission is optically thick the observed lu- 
minosity depends on the brightness temperature, Tb, as 
well as the line width, so that Leo on cjTb- Substituting 
in the relationship between density and line width. 



aco oc 



(12) 



The brightness temperature, Tb, will depend on the ex- 
citation of the gas (Eq. [6|) and the filling fraction of 
emission in the telescope beam, /;,. For high density and 
optical depth the excitation temperature will approach 
the kinetic temperature. Under those conditions, Eq. [T^ 
also implies aco « p°-^ifb Tkin)~^ ■ 

Thus, even for virialized GMCs we expect that the 
CO-to-H2 conversion factor will depend on environmen- 
tal parameters such as gas density and temperature. To 
some degree, these dependencies may offset each other. 
If denser clouds have higher star formation activity and 
are consequently warmer, the opposite effects of p and 
Tb in Eq. [12] may partially cancel yielding a conversion 
factor that is closer to a constant than we might other- 
wise expect. 

We also note that relation between mass and luminos- 
ity expressed by Eq. [TU] is not exactly linear, which is 
the reason for the weak dependence of aco on Leo or 
M^oi in Eq. [TTJ As a consequence, even for GMCs that 
obey this simple picture, aco will depend (weakly) on 
the mass of the cloud considered, varying by a factor of 
~ 4 over 3 orders of magnitude in cloud mass. 

2.2. Galaxies 

This simple picture for how the CO luminosity can be 
used to estimate masses of individual virialized clouds 
is not immediately applicable to entire galaxies. An ar- 
gument along similar lines, however, can be laid out to 
suggest that under certain conditions there should be 
an approximate proportionality between the integrated 
CO luminosity of entire galaxies and their molecular 
mass. This is known as the "mist" model, for rea- 
sons t hat will become clear in a few pa ragraphs. Fol- 
lowing iDickman. Snell k, Schloerbl (I1986D , the luminos- 



ity due to an ensemble of non-overlapping CO emit- 
ting clouds is Leo J^i^i'^Biai) cri, where is the 
area subtended by cloud i, and TBifli) and Oi are its 
brightness temperature and velocity dispersion, respec- 
tively. Under the assumption that the brightness tem- 
perature is mostly independent of cloud size, and that 
there is a well-defined mean, Tb, then Ts(ai) w Tg. 
We can rewrite the luminosity of the cloud ensemble 



as Leo ~ V2t:Tb Nciouds < TrUfaiRi) >, where the 
brackets indicate expectation value, Ndouds is the num- 
ber of clouds within the beam, and we have used = 
7ri?|. Similarly, the total mass of gas inside the beam is 
M„ioi « Ndouds < 4/37ri?fp(i?i) >, where p(i?^) is the 
volume density of a cloud of radius Ri. Using our defi- 
nition from Eq. [5] and dropping the i indices, it is then 
clear that 
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-'CO 



< R^pjR) > 

9n Tb < R^a{R) > ' 



(13) 



If the individual clouds are virialized they will follow 
Eq. [H or equivalently a = 0.0635 Ry/p. Substituting 
into Eq. [13] we find 



aeo oc 



< R^p > 



Tb < R^p" '' >' 



(14) 



which is an alogous to Eq. [T^ (obtained for i ndividual 
clouds). As IDickman. Snell fc Schloerbl ()1986[ ) discuss, 
it is possible to generalize this result if the clouds in a 
galaxy follow a size-line width relation and they have a 
known distribution of sizes. Assuming individually viri- 
alized clouds, and using the size-line width relation (Eq. 
[5]) , we can write down Eq. [T3] as 



aeo = 7.26 



CMC 



< i?2 > 



Te 



(15) 



where we have introduced the explicit dependence of the 
coefficient of the size-line width relation on the cloud 
surface density, Sqmc- This equation is the analogue of 
Eq. [n] 

In the context of these calculations, CO works as a 
molecular mass tracer in galaxies because its intensity is 
proportional to the number of clouds in the beam, and 
because through virial equilibrium the contribution from 
each cloud to the total luminosity is approximately pro- 
portional to its mass, as discussed for individual GMCs. 
This is the essence of the "mist" model: although each 
particle (cloud) is optically thick, the ensemble acts opti- 
cally thin as long as the number density of particles is low 
enough to avoid shadowing each other in spatial-spectral 
space. 

Besides the critical assumption of non-overlapping 
clouds, which could be violated in environments of very 
high density leading to "optical depth" problems that 
may render CO under luminous, the other key assump- 
tion in this model is the virialization of individual clouds, 
already discussed for GMCs in the Milky Way. The ap- 
plicability of a uniform value of aeo across galaxies relies 
on three assumptions that should be evident in Eq. 1151 a 
similar value for EgmCi similar brightness temperatures 
for the CO emitting gas, and a similar distribution of 
GMC sizes that determines the ratio of the expectation 
values < R^ > / < R^-^ >. Very little is known cur- 
rently on t he distribution of GMC sizes outside the Lo - 
cal Group (|Blitz et al.ll2007t iFukui k Kawamurall2010D , 
and although this is a potential source of uncertainty in 
practical terms this ratio is unlikely to be the dominant 
source of galaxy-to-galaxy variation in aeo • 

2.3. Other Sources of Velocity Dispersion 
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Because CO is optically thick, a crucial determinant 
of its luminosity is the velocity dispersion of the gas, 
a. In our discussion for individual GMCs and ensem- 
bles of GMCs in galaxies we have assumed that a is 
ultimately determined by the sizes (through the size- 
line width relation) and virial masses of the clouds. It 
is especially interesting to explore what happens when 
the velocity dispersion of the CO emission is related 
to an underlying mass distribution that includes other 
components besides molecular gas. Following the rea- 
soning by Downcs, Solomon & Radford (199 3,) (see also 
iMalonev fc Black! Il988l: iDownes fc Solo5ioi3ll99l and 
the discussion in ^2.11 we can write a cloud luminos- 
ity Leo for a fixed Tb as Lqq — Lqo o'*/(T, where the 
asterisk indicates quantities where the velocity disper- 
sion of the gas is increased by other mass components, 
such as stars. Assuming that both the molecular gas 
and the total velocity dispersion follow the virial ve- 
locity dispersion due to a uniform distribution of mass, 
a = ^jGM/bR, then = LQo\fM*jM^u where 

M* represents the total mass within radius R. Substi- 
t uting Lqo = MmoiloiQo yields the result aco-^co = 

Therefore, the straightforward application of aco to 
the observed luminosity Lqq will yield an overestimate 
of the molecular gas mass, which in this simple reason- 
ing will be the harmonic mean of the real molecular mass 
and the total enclosed mass. If the observed velocity dis- 
persion is more closely related to the circular velocity, as 
may be in the center of a galaxy, then a* « ^ GM* /R 
and the result of applying aco to Lqq will be an even 
larger overestimate of Mmoi- The appropriate value of 
the CO-to-H2 conversion factor to apply under these cir- 
cumstances in order to correctly estimate the molecular 
mass is 

where ^ is a geometrical correction factor accounting 
for the differences in the distributions of the gas and 
the total mass, so that K, = (a* ja)^ . In the extreme 
case of a uniform distribution of gas responding to the 
potential of a rotating disk of stars in a galaxy center, 
/C ^ 5. Everything else being equal, in a case where 
M* ^ 10 Af„io; the straight application of a standard 
aco in a galaxy center may lead to overestimating Mmoi 
by a factor of '-^ 7. 

Note that for this correction to apply the emission 
has to be optically thick throughout the medium. Oth- 
erwise any increase in line width is compensated by a 
decrease in brightness, keeping the luminosity constant. 
Thus this effect is only likely to manifest itself in regions 
that are already rich in molecular gas. Furthermore, it 
is possible to show that an ensemble of virialized clouds 
that experience cloud-cloud shadowing cannot explain a 
lower Xcoj simply because there is a maximum attain- 
able luminosity. Therefore we expect Xqo to drop in 
regions where the CO emission is extended throughout 
the medium, and not confined to collections of individual 
self-gravitating molecular clouds. This situation is likely 
present in ultra- luminous infrared galaxies (ULIRGs), 
where average gas volume densities are higher than the 



typical density of a CMC in the Milky Way, suggesting 
a per vading molecular ISM fe.g.. lScoville. Yun fc BrvantI 
I1997D . Indeed, the reduction of Xco in mergers and 
galaxy cente rs has been modeled in detail b v lShettv et al.l 
(|2011b ^ and'Naravanan et al.l (20111I2012D . and directly 
observed (gElEl SO gTJ) 

2.4. Optically Thin Limit 

Although commonly the emission from ^^CO 
J = 1 — > transition is optically thick, under con- 
ditions such as highly turbulent gas motions or 
otherwise large velocity dispersions (for example stellar 
outflows and perhaps also galaxy winds) then emission 
may turn optically thin. Thus it is valuable to consider 
the optically thin limit on the value of the CO-to-Il2 
conversion factor. Using Eq. [51 the definition of optically 
thin emission (Ij = tj [Bj{Tex) ~ Bj{Tcmb)], where 
Bj is the Planck function at the frequency vj of the 
J — !■ J — 1 transition, T^x is the excitation temperature, 
and Tcmb is the temperature of the Cosmic Microwave 
Background), and the definition of antenna temperature 
Tj, I J — {2ki'j/c^)Tj, the integrated intensity of the 
J — > J — 1 transition can be written as 

W(CO) = TjAv = ^^^,^±f,^,Nj. (17) 

The factor fcmb accounts for the effect of the Cos- 
mic Microwave Background on the measured intensity, 
/cmh = 1 - Z'"'^-- - l)/{e''''j/kT.^'' - 1). Note that 

fcmh ^ 1 fol ^ex '^cmb • 

The column density of H2 associated with this in- 
tegrated intensity is simply iV(H2) = SjLo -^J^' 
where Zqq is the CO abundance relative to molecu- 
lar hydrogen, Zqq = CO/II2. For a Milky Way gas 
phase carbon abundance, and assuming all gas-phase 
carbon is locked in CO molecules, Zqo ~ 3.2 x 10~^ 
([Sofia et al.ll2004[ ). Note, however, that what matters is 
the integrated Zqo along a line-of-sight, and CO may 
become optically thick well before this a bundance is 
reache d (for example. Fig. [T]). Indeed, Sheffe r et al.l 
([2008D analyze Zqo in Milky Way lines-of-sight, find- 
ing a steep Zco ~ 4.7 x IQ-\N {II2) / 1^'^^ cm-2)2.07 fo^. 
iV(Il2) > 2.5 X 10^° cm~^, with an order of magnitude 
scatter (see also ISonnentrucker et al.|[2007 f). 

When observations in only a couple of transitions 
are available, it is useful to assume local thermody- 
namic equilibrium applies (LTE) and the system is de- 
scribed by a Boltzmann distribution with a single tem- 
perature. In that case the column density will be 
iV(CO) Q(Te^)e'^i/'=^-Afi/gi, where Ei is the en- 
ergy of the J = 1 state [Ei/k « 5.53 K for CO), and 

Q{Tex) — ^2^=0 9 J corresponds to the parti- 
tion function at temperature T^x which can be approx- 
imated as Q{Tex) "-^ 2kTex/Ei for rotational transitions 
when Tex ^ 5.5 K ([Penziaslll975l . note this is accurate 
to ^ 10% even down to Tex ^ 8 K). Using Eq. [T7|we can 
then write 

_ iV(H2) ^^ 1 6h f Tex \ £iA 

^° W(CO) ^ Zco 8TT^fi^f„nbEi/k \Ei/kJ ■ 

(18) 
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Consequently, adopting Zco = 10 and using a repre- 
sentative Tex = 30 K, we obtain 

Xco « 1.6 X 10^^ ^ e^-°-i«4 cm-2(K km s-^)-\ 

(19) 

or aco ~ 0.34 M© (K km s ^ pc^) ^. These are an or- 
der of magnitude smaller than the typical values of Xqq 
and aco in the Milky Way disk, as we will discuss in 
21 Note that they are approximately linearly dependent 
on the assumed Zco and Tex (for Tex ^ 5.53 K). For 
a similar calc ulation that also includes a n expression for 
non- LTE, see lPapadopoulos et all (|2012D . 

2.5. Insights from Cloud Models 

A key ingredient in further understanding Xco in 
molecular clouds is the structure of molecular clouds 
themselves, which plays an important role in the radia- 
tive transfer. This is important both for the photodis- 
sociating and heating ultraviolet radiation, and for the 
emergent intensity of the optically thick CO lines. 

The CO J = 1 — >■ transition arises well within the 
photodissociation region (PDR) in clouds associated with 
massive star formation, or even illumi nated by the gen- 
eral d i ffuse interstellar radiation field (iMalonev fc Blac3 
[19881 IWolfire. Hollenbach fc Tielensl |1993D. At those 
depths gas heating is dominated by the grain photoelec- 
tric effect whereby stellar far-ultraviolet photons are ab- 
sorbed by dust grains and eject a hot electron into the 
gas. The main parameter governing grain photoelectric 
heating is the ratio x'^fcin/^e, where x is a measure of 
the far-ultraviolet field strength, and Ue is the electron 
density. This process will produce hotter gas and higher 
excitation in starburst galaxies. At the high densities of 
extreme starbursts, the gas temperature and CO exci- 
tation may also be enhanced by collisional coupling be- 
tween gas and warm dust grains. 

Early efforts to model the CO excitation and luminos- 
ity in molecular clouds usin g a large velo city gradient 
model were carried out by Goldre ich fc KwanI |.1974). 
Th e CO luminosity -gas mass relation was investigated 
by iKutner fc Leung ( 1985) using microtu r bulent mod- 
els, and~bjr^^IfireJJollenbach2^ ()1993D using 
both microturbulent and macroturbulent models. In mi- 
croturbulent models the gas has a (supersonic) isotropic 
turbulent velocity field with scales smaller than the pho- 
ton mean free path. In the macroturbulent case the scale 
size of the turbulence is much larger than the photon 
mean free path, and the emission arises from separate 
Doppler shifted emitting elements. Microturbulent mod- 
els produce a wide range of CO J = 1 profile shapes, 
including centrally peaked, flat topped, and severely cen- 
trally self-reversed, while most observed line profiles are 
centrally peaked. Macroturbulent models, on the other 
hand, only produce centrally peaked profiles if there are 
a sufficient number of "clumps" within the beam with 
densities n > 10'^ cm~'^ in order to pr o vide the peak 
brightness temperature. iFalgarone et all (|1994l ) demon- 
strated that a turbulent velocity field can produce both 
peaked and smooth line profiles, much closer to observa- 
ti ons than macroturbulent mode ls. 

IWolfire. Hollenbach fc Tieleni (fT993l) use PDR mod- 
els in which the chemistry and thermal balance was cal- 



culated self-consistently as a function of depth into the 
cloud. The microturbulent PDR models are very suc- 
cessful in matching and predicting the intensity of low- J 
CO lines and the e mission from many other atomic an d 
molecular species fH oUenbach fc Tielensl I1997L [1999"). 
For example, the nearly constant ratio of [CII]/C0 
J = 1 — >■ observed in both Galactic and extraga lactic 
sources ([Crawford et al.l fl98^ IStacev et al.l [T99ll) was 
first explained by PDR models as arising from high den- 
sity (n ^ 10"^ cm"'^) and high UV field (x ^ 10^) sources 
in which both the [CII] and CO are emitted from the 
same PDR regions in molecular cloud surfaces. These 
models show that the dependence on CO luminosity with 
incident radiation field is weak. This is because as the 
field increases the rco — 1 surface is driven deeper into 
the cloud where the dominant heating process, grain pho- 
toelectric heating, is weaker. Thus the dissociation of CO 
in higher fields regulates the temperature where rco = 1- 
PDR models have the advantage that they calculate 
the thermal and chemical structure in great detail, so 
that the gas temperature is determined where the CO 
line becomes optically thick. We note however, that al- 
though model gas temperatures for the CO J = 1 — )• 
line are consistent with observations the model tempera- 
tures are typically too c ool to match the ob served high-J 
CO line emission (;e.g.. lHabart et al.l[2010l ). The model 
density structure is generally simple (constant density or 
constant pressure), and the velocity is generally consid- 
ered to be a constant based on a single microturbulent 
velocity. More recent dynamical models have started to 
combine chemical and thermal calculations with full hy- 
drodynamic simulations. 

3. A MODERN THEORETICAL PERSPECTIVE 

In recent years there has been much progress in updat- 
ing PDR models, and in combining hydrodynamic simu- 
lations, chemical modeling, and radiation transfer codes. 

3.1. Photodissociation Regions 

Since CO prod uction mainly occ urs through ion- 
neutral chemistry ( van Dishoeck &: Bla ck 1988) the ion- 
ization structure through the PDR is important in set- 
ting the depth of the CO formation. Recombination of 
metal ions on PAHs can modify the abundance of free 
electrons. For a typical CMC with n ~ 10^ cm~^ il- 
luminated by a radiation field ~ 10 times the Galactic 
ISRF (the interstellar radiation field in the vicinity of the 
Sun) , the rco = 1 surface is at a depth of Ay ~ 1 when 
PAHs are included. Without PAHs the electron abun- 
dance stays high, and the ion-neutral chemistry is slowed 
due to Hj^ recombination. Consequently, the rco 1 
surface is pushed deeper in to the cloud (Ay ^ 2). The 
PAH rates are estimated bv IWolfire et al.l ([2008D by con- 
sidering the C°/C+ ratio in diffuse lines of sight , but 
there is considerable uncertainty both in the rates ([Liszti 
120 111 ) and in the PAH abundance and their variation with 
cloud depth. 

Bell et aQ ([2006( ) carry out a parameter study of Xco 
using PDR models with constant microturbulent line 
width. The calculated Xco versus Ay plots have a char- 
acteristic shape with high A"co at low Ay, dropping to 
a minimum value at Ay ~ 2 — 4 and then slowly ris- 
ing for increasing Ay. The Xco dependence at low Ay 
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arises from molecular gas with low CO abundance. Af- 
ter the CO line intensity becomes optically thick, Xqo 
slowly rises again as A^(H2) increases. Increasing density 
up to the critical density of CO J = 1 — ^ (ri(H2)cr,i ^ 
2200 cm~^) enhances the CO excitation and causes Xqo 
to drop. In addition, the minimum moves closer to the 
cloud surface as does the rco = 1 surface. An increase 
in Xco by a factor of more than 100 is found when de- 
creasing dust and metallicity to 1% of the local Galactic 
ISM, r eflecting the larger column of II2 at a given Ay. 
iBell. Viti fc Williams! (|2007,) suggest that the appropri- 
ate Xco value to use in various extragalactic environ- 
ments can be estimated from the minimum in the Xco 

v ersus Ay plot. 

iWolfire. Hollenbach fc McKei (20W) use an updated 
version of their PDR models, including a self-consistent 
calculation for the median density expected from a tur- 
bulent density distribution. The models provide a the- 
oretical basis for predicting the molecular mass fraction 
outside of the CO emitting region in terms of incident 
radiation field and gas metallicity. They assume that 
the median density is given by (n)mod = nexp(/i), where 
n is the volume averaged density distribution n oc 1/r, 
and u = 0.5 ln(l + 0.257W^) (jPadoan. Jones fc Nordlundl 
119971) . The sound speed that enters in the Mach number, 
Ai, is calculated from the PDR model output while the 
turbulent velocity is given by the size-linewidth relation 
(Eq.HD. 

3.2. Numerical Simulations 

As an alternative to PDR models with simple geome- 
tries and densities, hydrodynamical models can be used 
to calculate the line width, and density, but with only 
limited spatial res olution and approximate ch emistry and 
thermal balance. IGlover fc Mac Lowl (j2007l iH) carry out 
simulations of the formation of molecular clouds using 
a modified version of the magnetohydrodynamical code 
ZEUS-MP. They include a time-dependent chemistry, in 
particular for II2 formation and destruction, and ther- 
mal balance. IGlover et al.l (j2010 ') enhance the code to 
account for CO chemistry in a turbulent CMC. A turbu- 
lent picture of a CMC is not one with well defined clumps 
surrounded by an interclump medium but one with a 
continuous density distribution and constant mixing be- 
tween low and high density regions. The CO abundance 
varies within the cloud depending on the gas density and 
penetration of the dissociating radi ation field. Calcu- 
lation s of Xco are carri e d out by IGlover fc Mac Lowl 
(|20Tll) . and iShettv et al.) (|20Tll lH). The latter two ref- 
erences include non-LTE CO excitation and line trans- 
fer in the LVG approximation. These authors carry out 
3D turbulent simulations in a box of fixed size (20 pc), 
and turbulence generated with uniform power between 
wavenumbers 1 < fc < 2, but with various initial den- 
sities and metallicities. For their standard cloud model 
the saturation amplitude of the ID velocity dispersion 
is 2.4 km s^^. Each line of sight has a different CO in- 
tensity, and CO and H2 column density depending on 
the past and present physical conditions along it. Thus, 
for a given jV( H2) there is a rang e in Xco values (see 
Figs. 5 and 6 in IShettv et al.ll20Tl[ ). The dispersion gen- 
erally increases for clouds with lower density and lower 
metallicity. 



Figure 7 in IShettv et al.l (|2011[ ) shows the calculated 
mean Xco in different Ay bins for several initial densi- 
ties and metallicities. The variation in Xco with Ay is 
qualitatively similar to that shown in IBell et aD ()2006l ) 
for microturbulent models. Lower densities and metallic- 
ities drive Xco to higher values due to lower CO excita- 
tion and CO/II2 ratios respectively. At Ay > 7 the mod- 
els with different initial densities converge as the CO line 
becomes optically thick. The minimum in Xco occurs 
at larger Ay in the hydrodynamic simulations compared 
to the one sided PDR models, likely due to the dissoci- 
ating radiation incident on all sides of the box and due 
to its greater penetration along low density lines of sight. 
A range in density n = 100 — 1000 cm~'^ and metallicity 
Z = 0A — 1Zq produces a range of only Xco, 20 ~ 2—10. 
Thus although there can be large variations in Xco along 
different lines of sight in the cloud the emission weighted 
Xco is close to the typical Galactic value In ad- 

dition, the low metallicity case has higher Xco at low 
Ay but approaches the Galactic value for higher density 
(higher Ay) lines of sight. Although the models were run 
with constant box size, the dependence on Ay suggests 
that clouds with sufficiently small size so that CO does 
not become optically thick will have hig her Xqq . The re - 
sults seem to confirm the suggestion by IBell et al.l (|2006D 
that the mean Xco value should be near the minimum 
when plotted a s Xco v ersus Ay. 

She ttv et al.l (|2011bl ) investigate the results of vary- 
ing the temperature, CO abundance, and turbulent line 
width. They find only a weak dependence on temper- 
ature with Xco oc T-JJ-^ for 20 K< T^^ <100 K, and 
thus over the range of temperatures typically found in 
Galactic clouds Xco is not expected to vary signifi- 
cantly due to temperature. At low (fixed) CO abundance 
(nco/'^H2 10~^) , the line does not become optically 
thick and thus W(CO) follows X(H2) with constant Xco 
up to at least logX(H2) = 22.5. Varying the turbulent 
line width increases the CO line intensity as expected 
since decreasing the self-absorption allows more CO line 
emission to escape the cloud. They also find, however, a 
decreasing brightness temperature and over the range of 
velocity dispersion cr = 2 — 20 km s~^, W(CO) changes as 
W(CO) oc cr°'^ instead of changing linearly. The higher 
turbulent velocities create dense shocks but also larger 
voids of low density material. The competing effects 
produce a dependence on a that is slower than linear. 
Increasing the column density as well as the line width 
might produce a model result closer to W(CO) oc a. 

An interesting result is that Xco is not sensitive to 
the internal velocity profile within the cloud but only to 
the total line width, however that might be generated. A 
cloud need not obey a power law size-linewidth relation 
to have the same Xco as one with a Gaussian distribu- 
tion and identical dispersion velocity. Thus, clouds need 
not be "virialized" nor obey a size-linewidth relation to 
have the same Xco- The modeled Xco converges for 
high Ay clouds and thus the Xco factor is not expected 
to vary from cloud to clo ud as long as the there is a large 
enough column density. iShetty et al.l (|2011b( ) conclude 
that a nearly constant Xco is the result of the limited 
range in column densities, temperatures, and linewidths 
found in Galactic molecular clouds and that applying a 
constant Xco is approximately correct to within a factor 
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of ~ 2. 

Global models of Xco using hydrodynamic simulations 
to investigate the variation due to g alactic environment 
are carried out bvlFe ldmann. Gncdi n fc Kravtsovl (poTl ) 
and lNaravanan et alj (|2dliU2012[ ). A significant problem 
in numerical simulations is how to handle the physical 
conditions and/or line e mission within regions smaller 
than the spatial grid. iFeldmann. Gnedin fc Kravtsovl 
J2012) use high res o lution 0.1 pc) simulations from 
[Glover fc Mac Lowl (|2011| ) for "sub-grid" solutions to 
cosmological simulations of 60 pc resolution. These mod- 
els used constant Tkin — 10 K, LTE excitation for CO, 
with either constant CO line width or one proportional 
to Smoi"'^- By comparing results at their highest res- 
olutions with those at 1 kpc, and 4 kpc, they asses 
the effects of spatial averaging on Xco- The averag- 
ing tends to reduce the variation of Xco on A^(H2) and 
UV radiation field intensity. At greater than kilopar- 
sec scales, and H2 column densities between 10^^ cm^^ 
and 10^^ cm~^, Xco changes by only a factor of 2. 
They find essentially no variation in Xco with UV field 
strength betwe en 0.1 and 100 t imes the Galactic ISRF. 
|R;ldniann. Gn edin fc Kravtsovl i2012) do find a signif- 
icant variation with metallicity, with a scaling Xcooc 
^-0.5 Jqj. virial line widths. Note that t his relation is 
much shallower than found by, for example. iGenzel et al.l 
()2012l ji8.2p . The authors suggest that low metallicity 
high-redshift galaxies may not obey the same gas surface- 
density to star formati on relation observe d i n loca l disks. 

The results of iNarayanan et al.l (|2011l ) and 
iNaravanan et al.l (|20l3 i will be discussed in §7.21 
Here we note that their effective resolution of ~ 70 pc 
requires adopting sub-grid cloud properties. Although 
the surfaces of GMCs more massive than ~ 7 x 10^ 
Mq are resolved, their internal structure is not. The 
resulting line and continuum transfer, and temperature 
and chemical structure can only be approximate. The 
resolution problem also enters in simulations of individ- 
ual clouds. We showed in ^Jl]that CO becomes optically 
thick within a column iV(H2) « 2 - 3 x 10^° cm"2_ 
Thus for sub parsec resolution 0.1 pc) and density 
greater than n > 10'^ cm~^, the physical conditions 
(temperature, density, and abundances) are averaged 
over the line forming region and the calculated emitted 
intensity can be in error. The resolution problem is 
much more severe for H2 dissociation. Self-shielding of 
H2 starts within a column of only j V(H2) ^ 10^** cm~^ 
(|de Jong. Boland fc Dalgarnol I1980D . Thus with a 
resolution of ~ 0.1 pc, and density of n > 10 cm~'^, the 
optical depth to dissociating radiation is already r > 10^ 
in a resolution element. A complementary approach is 
to apply a state-of-the-art PDR code with well resolved 
H2 formation to the output of hydrodynamic simulations 
()Levrier et aLllMl . With increasing computing power 
the issues of resolution will continue to improve. We 
suggest that high priority should be placed on creating 
large scale galactic simulations that are well matched 
to small scale simulations with resolved cloud structure. 
The former provide environmental conditions and cloud 
boundary conditions while the latter provide the chem- 
istry and line emission in a realistic turbulent cloud. 
Expanding the library of GMC models with a range of 
column densities, line widths, and external heating, and 



thoroughly checking them against observations, would 
be most helpful. 

4. Xco IN THE MILKY WAY 

The Galaxy is the only source where it is possible to 
determine the CO-to-H2 conversion factor in a variety of 
ways. It thus provides the prime laboratory to investi- 
gate the calibration and the variations of the proportion- 
ality between CO emission and molecular mass. 

In the following sections we will discuss three types of 
Xco determinations: 1) employing virial masses, a tech- 
nique that requires the ability to spatially resolve molec- 
ular clouds to measure their sizes and kinematics, 2) tak- 
ing advantage of optically thin tracers of column density, 
such as dust or certain molecular and atomic lines, and 
3) using the diffuse 7-ray emission arising from the pion 
production process that takes place when cosmic rays 
interact with interstellar medium protons. Gamma-ray 
techniques are severely limited by sensitivity, and are 
only applicable to the Milky Way and the Magellanic 
Clouds. The good level of agreement between these ap- 
proaches in our own galaxy is the foundation of the use 
of the CO-to-Il2 conversion factor in other galaxies. 

4.1. Xco Based on Virial Techniques 

The application of the virial theorem to molecular 
clouds has been discu ssed by a number o f auth ors, and 
recently reviewed by iMcKee fc Ostriked ()2007[ ). Here 
we just briefly summarize the fundamental points. The 
virial theorem can be expressed in the Lagrangian (fixed 
mass) or Eulerian (fixed volume) forms, the latter par- 
ticularly applicable to turbulent clouds where mass is 
constantly exchanged with the surrounding medium. In 
the somewhat simpler Lagrangian form, the virial equi- 
librium equation is 

2{K - Ks) + B + W = (20) 

where K is the volume integral of the thermal plus ki- 
netic energy, Ks is the surface pressure term, B is the 
net magnetic energy including volume and surface terms 
(which cancel for a completely uniform magnetic field), 
and W is the net gravitational energy which is deter- 
mined by the self-generated gravitational potential if the 
acceleration due to mass external to the cloud can be ne- 
glected. In the simple case of a uniform, unmagnetized 
sphere virial equilibrium implies 2K + W = 0. It is useful 
to define the virial parameter, amr, which corresponds to 
the ratio of total kinetic energy to gravitational energy 
(|Bertoldi fc McKedlT99l . so that ayrr = 5Ra^/GM. 

4.1.1. Are Clouds Vmalized? 

In this context, gravitationally bound objects have 
ayir — 1 . Whether interstellar clouds are entities in virial 
equilibrium, e ven in a time or ensemble averaged sense 
(jMcKeel [19991 ). is a matter of current debate. Observa- 
tional evidence can be interpreted in terms of systems 
out of equilibrium with rapid star format ion and su bse- 
quent disruption in a few Myr (e.g., Elm egreenll2000l ) . an 
evolutionary progression and a typical lifetime of a few 
tens of Myr, long enough for clouds to become virialized 
(e.g.. lBlitz fc"Shulll980bl:lFukui fc Kaw amura 2010), or a 
lifetime of hundreds of Myr (e.g., lScoville fc Hershlll979l) . 
iRoman-Duval et al.l (|2010l) find a median a^ir « 0.5 for 
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clouds in the inner Galaxy, suggesting that they are 
bound entities where Myir represents a reasonable mea- 
sure of the molecular mass, although casti ng doubt on 
the as sumption of exact virial equilibrium. iWong et alJ 
estimate a very large scatter in a^ir in the 
Large Magellanic Cloud, but do lack an independent 
mass tracer so their results rest on the assumption of 
a fixed Xco ■ Observations in the outer Galaxy show an- 
other angle of the situation. iHever. Carpenter fc Snelll 
(|2001h find that clouds with M^m > 10^ Mq are self- 
gravitating, while small clouds with masses Mmoi < 10'^ 
Mq are overpressured with respect to their self-gravity, 
that is, have Umr 1 and are out of equilibrium. Given 
the observed mass function, however, such clouds repre- 
sent a very small fraction of the molecular mass of the 
Milky Way. 

In any case, observed GMC properties can be under- 
stood as a consequence of approximate energy equipar- 
tition, which observationally is very difficult to distin - 
guish from virial equilibrium ()Ballesteros-Paredesll2006l ) . 
Clouds with an excess of kinetic energy, ayir ^ 1, per- 
haps due to ongoing star formation or SNe would be 
rapidly dissipated, while clouds with a dearth of kinetic 
energy, <C 1, would collapse at the free-fall velocity 
which is within 40% of the equipartition velocity disper- 
sion and challenging to distinguish from turbulent mo- 
tions in observations. Furthermore, the resulting star 
formation will inject energy into the cloud acting to re- 
store the balance. Thus from the standpoint of determin- 
ing cloud masses over large samples, the assumption of 
virial equilibrium even if not strictly correct, is unlikely 
to be very wrong. 

4.1.2. Observational Results 

The most significant study of the relation between 
virial mass and Lqo (the mass- luminosi t y rela tion) in 
the Milky Way is that bv iSolomon et all (|1987| ). which 
encompasses 273 clouds and spans several orders of mag- 
nitude in cloud luminosity and mass. It is dominated by 
clouds located in the inner Galaxy, in the region of the 
so-called Molecular Ring, a feature in the molecular sur- 
face density of the Milky Way peaking at Roc ~ 4 kpc 
galactocentric radius. It uses kinematic distances with 
an old value of the distance to the Galactic Center, 
Rq = 10 kpc. We report new fits after a 0.85 scal- 
ing in all distances and sizes and 0.72 in luminosities 
to bring them into agreement with the modern distance 
scale {Rq = 8.5 kpc). The virial mass computations 
a ssume a p(r) oc (s ee ^2.ip . 

ISolomon et al.l (|1987[ ) find a very strong correlation be- 
tween Myir and Leo, such that Myir = 37.9 Lco°'*^ with 
a typical dispersion of 0.11 dex for Myir. Note the excel- 
lent agreement with the expected mass-luminosity rela- 
tion in Eq. [lOJ using a typ ical CO brightness temperature 
Ts 4 K ((Malonevll 19901) . For a cloud at their approxi- 
mate median luminosity, Leo ~ 10^ K km s~^ pc^, this 
yields aco = 4.6 Mq (K km s~-^ pc^)~^ and Xco,20 = 
2.1. Because the relation is not strictly linear aco will 
change by ^ 60% for an order of magnitude change in 
luminosity (Fig. [5]). Therefore GMCs with lower lumi- 
nosities (and masses) will have somewhat larger mass-to- 
light ratios and conversion factors than more luminous 
GMCs. 
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Fig. 2. — Relation between virial aco and CO luminosity for 
GMCs in the Milky Way (Solomon et al. 1987^). We have corrected 
the numbers in the original table to reflect the updated distance to 
the Galactic Center of 8.5 kpc. The dependence of aco ori Leo 
arises from the fact that the correlation between M^ir and Lqq has 
a nonlinear slope (Mvir Lqq'-^'^^^^'^'^^^), following the expec- 
tations from Eq. 1101 for approximately constant brightness tem- 
perature. This results in aco ~ 4.61 (Lco/lO^)"" , denoted 
by the red thick dashed line (the dispersion around this relation 
is ±0.15 dex). The nominal value at Lco = 10^ K km s^^ pc^ is 
illustrated by the thin black dashed line. 

Independen t analy sis using the same survey by 
IScoville et al.l (|1987D yields a very similar mass- 
luminosity relation. After accounting for the differ- 
ent coefficients used for the calculation of the virial 
mass, the relation is Myir — 33.5 ico°'^^- For a 
^co ~ 10^ K km s^^ pc^ cloud this yields aco — 
6.0 Mq (K km s"^ pc^)"! and Xco,20 = 2.8 (this work 
uses Rq — 8.5 kpc). Interestingly, there is no substantial 
difference in the mass-luminosity relation for GMCs with 
or without HII regions (Scovillc & Good 1989), although 
the latter tend to be smaller and lower mass, and have 
on average half of the velocity-integrated CO brightness 
of their strongly star-forming counterparts. The result- 
ing difference in Tg could have led to a displacement in 
the relation, according to the simple reasoning leading to 
Eq. 1101 but it appears not to be significant. 

4.1.3. Considerations and Limitations 

Besides the already discussed applicability of the virial 
theorem, there are a number of limitations to virial stud- 
ies. Some are practical, while others are fundamental to 
the virial technique. On the practical side, virial studies 
are sensitive to cloud definitions and biases induced by 
signal-to-noise. These will impact both the values of R 
and a used to compute the mass. In noise- free measure- 
ments isolated cloud boundaries would be defined using 
contours of zero emission, when in reality it is neces- 
sary to define them u sing a higher contour (for example, 
ISolomon et"al] 119871 u se a T g ~ 4 K CO brightness 
contour). IScoville et al.l (|1987D discuss the impact of this 
correction, studying the "curve of growth" for R and a 
as the definition contour is changed in high signal-to- 
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noise observations. Moreover, isolated clouds are rare 
and it is commonly necessary to disentangle many par- 
tially blended features along the line of sight. To measure 
a size clouds need to be resolved, and if appropriate the 
telescope beam size needs to be deconvolved to estab- 
lish the intrinsic cloud size. This is a major concern in 
extragalactic studies, but even Galactic datasets are fre- 
quently undersampled which affects the reliability of the 
R and Leo determinations. Given these considerations, 
it is encouraging that two comprehensive studies using 
independent analysis of the same survey come to values 
of aco that differ by only ^ 30% for clouds of the same 
luminosity. 

A fundamental limitation of the virial technique is that 
CO needs to accurately sample the full potential and size 
of the cloud. For example, if because of photodissocia- 
tion or other chemistry CO is either weak or absent from 
certain regions, its velocity dispersion may not accurately 
reflect the mass of the cloud. This is a particular con- 
cern for virial measurements in low metallicity regions 
(see although most likely it is not a limitation in 
the aforementioned determinations of Xqq in the inner 
Galaxy. 

4.2. Column Density Determinations Using Dust and 
Optically Thin Lines 

Perhaps the most direct approach to determining the 
H2 column density is to employ an optically thin tracer. 
This tracer can be a transition of a rare CO iso topologue 
or other chemical species (e.g., CH lMagnanlet al. 2003). 
It can also be dust, usually optically thin in emission at 
far-infrared wavelengths, and used in absorption through 
stellar extinction studies. 

4.2.1. CO Isotopologues 

A commonly used isotopologue is ^^CO. Its abun- 
dance relative to ^^CO is down by a factor approach- 
ing the ^^C/^'^C~ 69 isotopic ratio at the solar circle 
(^^C/^'^C« 50 at Rcc ~ 4 kpc, the galactocentric radius 
of the Molecular Ring) as long as chemical fractionation 
and selective photodissociation effects can be neglected 
(pVilson 1999). Given this abundance ratio and under the 
conditions in a dark molecular cloud ^^CO emission may 
not always be optically thin, as ti ~ 1 requires Ay ~ 5. 

The procedure consists of inverting the observed inten- 
sity of the optically thin tracer to obtain its column (or 
surface) density. In the case of isotopologues, this column 
density is converted to the density of CO using the (ap- 
proximate) isotopic ratio. Inverting the observed inten- 
sity requires knowing the density and temperature struc- 
ture along the line of sight, which is a difficult problem. If 
many rotational transitions of the same isotopologue are 
observed, it is possible to model the line of sight column 
density using a number of density and temperature com- 
ponents. In practice an approximation commonly used 
is local thermodynamic equilibrium (LTE), the assump- 
tion that a single excitation temperature describes the 
population distribution among the possible levels along 
the line of sight. It is also frequently assumed that ^^CO 
and ^'^CO share the same T^x, which is particularly jus- 
tifiable if collisions dominate the excitation {T^x = Tkin, 
the kinetic temperature of the gas) . Commonly used ex- 
pressions for determining iV(^^CO) under these assump- 
tions can be found in, for example, iPineda et al.. (,2010, ). 



Note, however, that if radiative trapping plays an im- 
portant role in the excitation of ^^CO, Tf-x for ^'^CO will 
generally be lower due to its reduced optical depth (e.g., 
IScoville fc Sandersl [l987bh. 

'Dic kmanI ()1978[ ) characterized the CO column density 
in over 100 lines of sight toward 38 dark clouds, focus- 
ing on regions where the LTE assumption is unlikely 
to introduce large errors. The combination of LTE col- 
umn densities with estimates of Ay performed using star 
counts yields Ay « (4.0±2.0) x IQ-^^N{^^C0) cm^ mag. 
Comparable re sults were obtained in det a iled s tudies 
of Taurus by iFrerking. Langer fc WilsonI (|1982l . note 
the nonlinearity in their express ion) and Perseus by 
IPineda. CaseUi 'fc GoodmanI (|2008D . t he latter usmg a so- 
phist icated extinction determination (jLombardi fc Alvei 
1200 It ). Extinction can be converted into molecular col- 
umn density, thr ough the assumption of an e ffective 
gas-to-dust ratio. iBohlin. Savage fc Drake! (jl97 8l) deter- 
mined a relation between column density and reddening 
(selective extinction) such that [N(m)+2N{Yi2)\/ E{B- 
V) ~ 5.8 X 10^^ atoms cm~^ mag""'^ in a survey of in- 
terstellar Lya absorption carried out using the Coperni- 
cus satellite toward 75 lines of sight, mostly dominated 
by HI. For a "standard" Galactic interstellar extinction 
curve with Ry = Ay /E{B — V) = 3.1, this results in 

Nh = N{m) + 2 iV(H2) w 1.9 X 10^^ cm-2 Ay. (21) 

A much more recent study using Far Ultraviolet Spectro- 
scopic Explor er observations fi nds essentially the same 
relation (Rac hford et al.|[2009[ ). In high surface density 
mole cular gas Ry may be closer to 5.5 ([Chapman et al.l 
l2009ll. and Eg. M may yield a 40% overestimate 
(jEvans et all |2009). Using Eq. EH the approxi- 
mate relation between ^^CO J = 1 — >■ and molecu- 
lar column density is AT (Ho) » 3.8 x 10^ N(}^CO). 
IPineda. Caselli &: GoodmanI (|2008 ^ find a similar result 
in a detailed s tudy of Perseus, with an increased scatter 
for Ay > 5. IGoldsmith erall f2008) use these results 
together with an averaging method to increase the dy- 
namic range of their ^'^CO and ^^CO data, a physically 
motivated variable ^^CO/^"^CO ratio, and a large veloc- 
ity gradient excitation analysis, to determine H2 column 
densities in Taurus. They find that Aco,20 ~ 1-8 recov- 
ers the molecular mass over the entire region mapped, 
while there is a marked increase in the region of low 
column density, where Xqo increases by a factor of 
5 where A^(H2) < 10^^ cm~^. As a cautionary note 
about the blind use of ^'^CO LTE estimates, however, 
iHeiderman et al.] (j2010( ) find that this relation between 
H2 and ""^^CO underestimates A^(H2) by factors of 4 — 5 
compared with extinction-based results in the Perseus 
and Ophiuchus molecular clouds. 

4.2.2. Extinction Mapping 

Extinction mapping by itself can be directly employed 
to determine Xqq- It fundamentally relies on the as- 
sumption of spatially uniform extinction properties for 
the bands employed, and on the applicability of Eq. [211 
t o convert extinction into col umn density. 

IFrerking. Langer fc: WilsonI (|1982D determined 
-^^00,20 ~ 1-8 in the range 4 < Ay < 12 in p Oph, while 
the same authors found constant W(CO) for Ay > 2 
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in Taurus. iLombardi. Alves fc Ladal (|2006[ ) studied the 
Pipe Nebula and found a best fit Xqo in the range 
-^CO,20 ~ 2.9 — 4.2, but only for K-band extinctions 
Ak > 0.2 (equivalent to Ay > 1.8, iRieke fc Lebofskvl 
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|1985() . A simple fit to the data ignoring this nonlinearit 
ields Xco,20 2.5. The Pine da. Caselli fc Goodma: 
study of Perseus finds Xco,20 ~ 0.9 — 3 over a 
number of regions. The relation between CO and H2, 
however, is most linear for Ay < 4, becoming saturated 
at larger line-of-sight extinctions. 




E 
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Fig. 3. — Relation between C O column density a nd extinction 
in the Taurus molecular cloud llPineda et al.1 [20101 ). The figure 
shows the pixel-by-pixel relation between gas-phase CO column 
density (obtained from '^^CO) and Ay- The blue line illustrates 
the "average" linear relation for 3 < Ay < 10, N{^-^CO) ~ 
1.01 X 10^"^ Ay cm-2 (implying CO/H2 Ri 1.1 X 10"* for the as- 
sumed isotopic ratio ). The linearity is clearly broken for Ay > 10. 
llPineda et al.|[2010l) show that linearity is restored to high Ay af- 
ter applying a correction for CO freeze-out into dust grain mantles. 



( 2010 ) extend the aforementioned 
[Goldsmith et al.l (2008) study of Taurus by character- 
izing the relation between reddening (from the Two 
Micron All Sky Survey, 2MASS) and CO column density 
(derived from ^'^CO) to measure Xco,20 ~ 2.1. They 
find that the relation between Ay and CO flattens 
for Ay > 10 (Fig. [3]), a fact that they attribute to 
freeze-out of CO onto dust grains causing the formation 
of CO and CO2 ice mantles. Including a correction for 
this effect results in a linear relation to Ay < 23. For 
Ay < 3 the column density of CO falls below the linear 
relationship, likely due to the effects of photodissoci- 
ation and chemical frac tionation. Along similar lines, 
iHeiderman et al.l ()2010( ) find that in Ophiuchus and 
Perseus CO can underpredict H2 with respect to Ay for 
E moi > 200 MfT) pc~^ by as much as ~ 30%. 

iParadis et al.l ()2012[ ) recently used a high-latitude ex- 
tinction map derived from 2MASS data using an exten- 
sion of the NICER methodology (Dobashi et al. 20Q1 
|2009( ) to derive Xco in sample of nearby clouds with 
l&l > 10°. They find Xco,20 » 1-67 ± 0.08 with a some- 
what higher value Xco, 20 ~ 2.28 ± 0.11 for the inner 
Galaxy region where |^| < 70°. They report an excess in 
extinction over the linear correlation between total gas 



and Ay at 0.2 < Ay < 1.5, suggestive of a gas phase 
that is not well traced by either 21 cm or CO emission. 
We will return to this in ^14.2.41 

4.2.3. Dust Emission 

The use of extinction mapping to study A^(H2) is 
mostly limited to nearby Galactic clouds, since it needs a 
background stellar distribution, minimal foreground con- 
fusion, and the ability to resolve individual stars to deter- 
mine their reddening. Most interestingly, the far-infrared 
emission from dust can also be employed to map the gas 
distribution. Indeed, dust is an extraordinarily egalitar- 
ian acceptor of UV and optical photons, indiscriminately 
processing them and reemitting in the far-infrared. In 
principle, the dust spectral energy distribution can be 
modeled to obtain its optical depth, Td{X), which should 
be proportional to the total gas column density under 
the assumption of approximately constant dust emissiv- 
ity per gas nucleon, fundamentally the product of the 
gas-to-dust ratio and dust optical properties. 

How valid is this assumption? An analysis of the corre- 
lation between and H I was carrie d out at high Galactic 
latitudes bv iBoulanger et al.l (|1996[ ). who found a typical 
dust emissivity per H nucleon of 
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Td/Nn w 1.0 X 10"^^(A/250Aim) 



-^cm2, 



(22) 



with /3 = 2, in excellent accord with the recent value 
for h igh latitude gas derived using P lanck observa- 
tions (jPlanck Collaboration et al.ll2011(l who prefer /3 = 
1.8). They also identified a break in the correlation for 
iV(HI) ^ 5 X 10^" cm~^ suggestive of an increasingly im- 
portant contribution fro m H2 to Nh. in agre ement with 
results from Copernicus (| Savage et al.l [19771 ). There is 
evidence that the coefficient in Eq. [22| changes in molec- 
ular gas. It may i ncrease by factors of 2 — 3 at very high 
column densities fSchnce et al."2008'; 'Fla gev et aLll2009l: 
Planc k Col laboration et al. 2011d), likely due to grain 
grow th or perhaps sol id state effects at low temperatures 
(e.g.. lMenv et al.ll2007f ). Note, however, that recent work 
using Planck in the Galactic plane finds Sbgr ~ (0.92 ± 
0.05) X 10~^^ cm^ at 250 ^ m, with no significant vari- 
ation with Galactic radius (iPlanck C ollaboration et al.l 
l2011b( ). This SuGR is almost identical to that observed 
in dust mixed with mostly atomic gas at high latitudes, 
suggesting that the aforementioned emissivity variations 
are very localized. 

This excellent correlation between and Nh is the 
basis for a number of studies that use dust emis- 
sion t o determine H2 column dens i ties. Most no- 
tably, iDame. Hartmann fc ThaddeusI (|2001[ ) employed 
the Columbia survey of molecular gas in the Galactic 
plane together with the Dwingeloo-Leiden HI survey and 
the IRAS temperature-corrected 100 /im spectral density 
map by Schlcgel, Finkbcincr & Davis (1998). With these 
data, iV(H2) can be obtained using 



iV(H2) - [Td/SuGK - A^(HI))/2, 



(23) 



which simply states that the dust optical depth (r^) is 
a perfect tracer of the total column density of gas when 
the emissivity per nucleon ((5dgr,) is known. As we just 
discussed, i^dgr can be straightforwardly determined on 
lines of sight dominated by atomic gas, for example. The 
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Fig. 4. — Planck results in the Aquila-Ophiuchus flare 
UPlanck Collaboration et al.|[2011al) . The figure shows molecular 
gas column density with a color range A'^(H2) ~ —1.5 X 10^^ 
cm~^(dark blue) to 3.5 X lO^'^ cm~'^(red). The black contours 
illustrate the CO emission from the Columbia survey, with values 
Ico ^ 2, 10, 20 K km s~^. The black region lacks CO information. 

comparison of the molecular column density so derived 
with the observed W(CO) yields Xco,20 ~ 1.8±0.3, vahd 
for |6| > 5° on large scales across the Galaxy. This study 
also finds evidence for a systematic increase in Xqo by 
factors ~ 2 — 3 at high Galactic latitude {b > 20°), in 
regions with typical iV(H2) < 0.5 x 10^° cm"^ on 0.5° 
angular scales. 

The excellent Planck dataset has afforded a new view 
on this topic (Fig. |4]). Planck C ollabo ration et al.l 
(|2011af) produced a new map for the Milky Way and 
a new determination of Xqo using Eq. [23] for lines of 
sight with l&l > 10°. They obtain Xco,20 = 2.54 ±0.13, 
somewhat larger than the previous study. This difference 
is likely methodological, for example the use of different 
far-infrared wavelengths as well as local versus global cal- 
ibrations of (5dgr- 

4.2.4. CO-Famt Molecular Gas and Diffuse Lines of Sight 

Most interestingly, IPlanck Collaboration et all (j2011af) 
observe a tight linear correlation between Td and A^(HI) -|- 
2XcoW(CO) for Av < 0.4 and 2.5 < Ay < 10, with 
an excess in in the intermediate range (Fig. [5]). 
This excess can be understood in terms of a compo- 
nent of H2 (or possibly a combination of H2 and cold, 
opaque HI) that emits weakly in CO and is prevalent at 
0.4 < Ay ^ 2.5 (the explanation is not unique, since 
the methodology cannot distinguish it from a change 
in dust emissivity in that narrow Ay regime, but such 
possibility appears unlikely). This molecular component 
arises from the region in cloud surfaces where gas is pre- 
dominantly H2 but most carbon is not in CO molecules 
because the extinction is too low, essentially the FDR 
surface (see ^6.11 and Fig. [TJ . This component is fre- 
quently referred to as "C O-dark molecular gas" or some- 
times simply "d ark gas" (^Grenier. Casan diian fc Terrieri 
[200I IWolfire. Hollenbach fc M cKcc 2010). In this re- 
view we will refer to it as "CO-faint," which is a more 
accurately descriptive name. The existence of molecular 
gas with low CO abundan ce has been noted previousl y 
in theoretical models (e.g.. Ivan Dishoeck fc BlacgiTOSSl) . 
in observations of diffuse gas and high-latitude clouds 



(e.g.. lLada fc Blit3ll 988D. and in obse rvations of irregu- 
lar galaxies ( e.g., [M adden et al. 19971). In this c ontext, 
the results by IPlanck Collaboration et al.l (|2011aD are in 
qualitative agreement with already discussed observa- 
tions that show an increase in Xgo at low molecular col- 
umn densities fe.g.. iGoldsmith et al.ll2008l : IParadis et al.l 
I2OII) . 




Fig. 5. — Correlation between t^j at 350 /im and total hy- 
drogen column density for Xqq 20 = 2.3 (for |6| > 10°, 
Planck Collaboration et al.|[2011al ). The color scale represents the 
logarithm of the number of lines-of-sight, and the blue dots the re- 
sult of Nil binning. The red dashed lines indicate Ay = 0.37 and 
Ay = 2.5. The red solid line represents the best linear fit for low 
A^H. Note that it is also a good fit to the Ay > 2.5 points. The ex- 
cess in the binned correlation over the red line for 0.37 < Ay < 2.5 
is either an indication of "CO-faint" molecular gas, or possibly a 
combination of high optical-depth "opaque" Hi with "CO-faint" 
H2, or a change in the dust emissivity over that Ay regime. 

The Planck observatio ns are also in qualitative agree- 
ment with the analysis bv lGrenier. Casandiian fc Terrieii 
(|2005f ). who correlated the diffuse 7-ray emission over 
the entire sky with templates derived from the HI, CO, 
and dust, finding a component of gas not traced by 
CO evident in local clouds at high latitudes. Their 
analysis finds that this component is as important, by 
mass, as the "CO-bright" H2 component in several of 
these clouds, and increasingly more important for smaller 
clo ud masses. The rec ent a nalyses based on Fermi dsda 
bv lAbdo et all (|2010df ) and I Ackermann et all (|2012f ) are 
also qualitatively compatible with these results, finding 
that the "CO-faint" component amounts to 40% - 400% 
of the "CO-bright" mass in the Cepheus, Polaris, Chama- 
leon, R Cr A, and Cassiopeia clouds (small local molec- 
ular clouds). 

The ionized carbon far-infrared fine structure emission 
provides an additional probe of molecular gas at low Ay ■ 
Large scale [CII] observations of the (^P3/2 — i'^Pi/2) fi^e- 
structure transition in the Milky Way and external galax- 
ies suggest it is due to a combination of emission from the 
Cold Neutral Medium ( CNM) and from PDRs located in 
the surfac es of GMCs (IStacev et all 119911 : IShibai et al] 
[1991; Be nnett et al.lll994D . The contribution from [CII] 
in the diffuse ionized gas, however could also be im- 
portant (|Heileslll99i [Madden et all I1993D parti cularly 
along certain lines-of-sight ([Velusamv eTaLllMll . 

In regions where most of the emission arises in PDRs 
[CII] has the potential to trace the "CO-faint" molec- 
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ular regime at low Ay- ILanger et alJ (|2010D analyze 
16 lines-of-sight in the plane of the Galaxy and find 
that in about half of them the observed [CII] inten- 
sity can be entirely explained as due to carbon in 
atomic gas in the CNM. The other half, however, ex- 
hibits [ClI]/iV(HI) ratios that are too large to be due 
to atomic gas and may have molecular to atomic ra- 
tios as large as iV(H2)/A^(Hl) - 6. While the very 
brightest [CII] components investigated arise from dense 
(n > 10^ cm~^) PDRs exposed to intense radiation fields, 
most of the [CII] emission in these molecular lines-of- 
sight can be explained as originating in the surfaces of 
modestly dense GMCs (n [3 - 300] x 10^ cm'^) ex- 
posed to at most a few times th e local interstellar ra - 
diation field at the solar circle (Pineda et al.' '2 010b( ). 
Fol lowing the reasoning in Vclusamy et al. (201(|) and 
ILanger et al.l (|2010( ). a very approximate relation be- 
tween [CII] emission and "CO-faint" H2 is A^(H2) ^ 
1.46 X 1020 W( [CII]) - 0.35A^(HI) cm'^, for W([CII]) in 
K km s~'^ (see Eq. [7] to convert between Jy and K). 
We caution that the coefficients correspond to the Milky 
Way carbon abundance and are very dependent on the 
assumed physical conditions, particularly the densities 
(we use n(HI) ~ 200 cm"^, ^(Ha) 300 cm'^, and 
Tht-Th2 - 100 K). 

ILiszt. Petv fc LucasI (poTot ) measure Xco in diffuse 
gas by first estimating the total hydrogen column 
from dust continuum emission (using the map by 
iSchlegel. Finkbeiner fc Davislll998D , and subtracting the 
observed HI column density (c.f., Eq. [23l). They se- 
lect lines-of-sight with HCO+ absorption spectra against 
bright extragalactic continuum sources. The CO emis- 
sion, together with the measured H2 column provides 
a mea sure of Xgo- Su r prisin gly, iLiszt. Petv fc LucasI 
(|2010D and ILiszt fc Petvl ()2012D find mean values in dif- 
fuse gas similar to those in GMCs. There are large vari- 
ations, however, about the mean with low Xco (bright 
CO) produced in warm T ^ 100 K diffuse gas and high 
Xco (faint CO) produced at low N(H2) column densi- 
ties. The authors argue that these variations mainly re- 
flect the CO chemistry and its dependence on ultraviolet 
radiation field, density, and total column density, rather 
than the H2 column density. We note that the observed 
CO column densities c annot be produced in steady- 
state PDR models (.e.g.. lSonnentrucker et aLi r2007l. En- 
hanced CO production might occur through the "CH+ 
channel" driven by non-thermal ion- ne utral reactions 
(jFederman et al.l 119961: iVisser et al']|2009D or by pockets 
of warm gas and ion-neutral reactions in turbulent dissi- 
pation regions (Godard, Falgarone fc Pincau dcs ForeJi 
|2009() . Density fluctuations in a tur bulent median migh t 
also increase the CO production (Levrier et al.l [20121 ). 
Thus, although the mean Xco diffuse cloud value is simi- 
lar to GMCs, the CO emission from diffuse gas cannot be 
easily interpreted as a measure of the molecular column 
except perhaps in a statistical sense. 

4.3. Xco Based on Gamma-Ray Observations 

Diffuse 7-ray emission in the Galaxy is chiefly due to 
three processes: neutral pion production and subsequent 
decay in collisions between cosmic-rays and interstellar 
matter, bremsstrahlung emission due to scattering of 
cosmic-ray electrons by interstellar matter, and inverse 



Compton scattering of low energy photons by cosmic- 
ray electrons. The first of these processes is the domi- 
nant production channel for diffuse 7-rays with energies 
above 200 McV, although at high Galactic latitude there 
will be an increasingly important inverse Compton com- 
ponent (iBloemen 1989). Interestingly, the fact that in- 
teractions between cosmic rays and nucleons give rise to 
diffuse 7-ray emission can be used to count nucleons in 
the ISM, and indeed the use of Xco to represent the ra- 
tio A^(H2)/W(C0) was introduced for the first time in 
-y-ra y work usi ng observations from the COS B satellite 
iLebrun et al.l 



Accounting for the pion decay and bremsstrahlung pro- 
cesses, and neglecting the contribution from ionized gas, 
the basic idea behind modeling the emission is to use 
a relation = J2 e^MR^) [N{m), + 2 Yco W(CO)J 
(jBloemenI 119891 ). Here I-y is the diffuse 7-ray emis- 
sion along a line-of-sight, is the HI gas emissiv- 
ity (a function of Galactocentric radius R), and Yco 
is a parameter that takes into account that emis- 
sivity in molecular clouds may be different than in 
atomic gas, due to cosmic-ray exclusion or concentration, 
Xco = Yco i-y .m / e-y m (.Gabici, A haronian fc Blasj 
l2007t iPadovani. Galli fc Glassgoldll20lj9[ ). This situation 
is analogous to that presented in dust emission tech- 
niques, where emissivity changes in the molecular and 
atomic components will be subsumed in the resulting 
value of Xco ■ 

Recent analyses use the approach qHiA^(HI) + 

gcoW(CO) -I- qEBvE{B - V)res + ■ ■ ■ , whcrc the first 
three terms account for emission that is proportional 
to the HI, the "CO-bright" H2, and a "CO-faint" H2 
component that is traced by dust reddening or emis- 
sion residuals (the additional terms not included ac- 
count for an isotropic 7-ray backg round and the con - 
tribution from point sources, e.g.. lAbdo et al.ll2010d[ ). 
The dust residual template consist s of a dust map (e.g., 
ISchlegel. Finkbeiner fc Davislll998l ) with a linear combi- 
nation of A^(HI) and W(CO) fitted and removed. 

4.3.1. Observational Results 

A disc ussion of the re sults of older analyses can be 
found in iBloemenI ()1989() , here we will refer to a few of 
the more recent results using the Compton Gamma Ray 
Observatory and Fermi sa tellites. 

iStrong fc Mattoxl (|1996( ) analyzed the EGRET all-sky 
survey obtai ning Xno.2n ^ 1.9±0 .2. Similar results were 
obtaine d bv iHunter et al.l (Il997f) for t he inn er Galaxy, 
and by iGrenier. Casandiian fc Terrieii ( 20051) for clouds 
in the solar neighborhood (see also >j4.2 4p . [Strong et al.l 
(|2004l) introduce in the analysis a Galactic gradient in 
Xco I in an attempt to explain the discrepancy between 
the derived e-y (R) and the distributions of supernova rem- 
nants and pulsars, which trace the likely source of cos- 
mic rays in supernovae shocks. Matching the emissivity 
to the pulsar distribution requires a significant change in 
Xco between the inner and outer Galaxy (Xco 20 ~ 0.4 
for R - 2.5 kpc to Xco,20 ~ 10 for i? > 10 kpcj. 

The sensitivity of Fermi has been a bo on for stud- 
ies of d iffuse 7-ray emission in our Galaxv. lAbdo et al.l 
l! 2010d[ ) and lAckermann eFall (|20Tll |2012b|(dD analyze 



the emission in the solar neighboorhood and the outer 
Galaxy. Taken together, they find a similar Xco for 
the Local arm and the interarm region extending out to 
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R ~ 12.5 kpc, Xco,20 ^ 1-6 — 2.1 depending on the as- 
sumed HI spin temperature. The very local high-latitude 
clouds in the Gould Belt have a lower J'^co,20 ~ 0.9, and 
appear no t to represent the averag e properties at the so- 
lar circle (jAckermann et aLll2012 b'f?). 

It is important to note, however, that the Fermi 
studies follow the convention Xco = Qco/'^Qm with 
the parameters defined in the previous section. Thus, 
their definition of Xqo docs not include the "CO- 
faint" envelope that is traced by the dust residual tem- 
plate: it only accounts for the H2 that is emitting 
in CO. This is different from the convention adopted 
in the dust studies discussed in the previous section, 
where all the H2 along a line-of-sight is associat ed with 
the corresponding CO. lAckermann et all (|2012b[) report 
^^1^00,20 « 0.96, 0.99, 0.63 for Chameleon, R CrA, and 
the Cepheus/Polaris Flare region, respectively using the 
Fermi convention. The ratio of masses associated with 
the CO and the dust residual template (reported in their 
Table 4) suggests that we should correct these numbers 
by factors of approximately 5, 2, and 1.4 respectively to 
compare them with the dust modeling. This correction 
results in Xco, 20 ~ 4.8, 2, 0.9, bringing these clouds in 
much better (albeit not complete) agreement with the 
results of dust modeling and the expectation that they 
be underluminous — not overluminous — in CO due to 
the presence of a large "CO-faint" molecular component. 

4.3.2. Considerations and Limitations 

We have already mentioned a limitation of 7-ray stud- 
ies of Xco , the degree to which they may be affected by 
the rejection (or generation) of cosmic rays in molecular 
clouds. A major limitation of 7-ray determinations is 
also the poor angular resolution of the observations. An- 
other source of uncertainty has become increasingly ap- 
parent with Fermi, which resolves the Magellanic Clouds 
m 7-ray emission (|Abdo et alllMoUbl ). In the Clouds, 
the distribution of emission does not follow the distribu- 
tion of gas. Indeed, the emission is dominated by regions 
that are not peaks in the gas distribution, but may corre- 
spond to sites of cosmic ray injection (e.g.. lMurphv et aLl 
[2012^ . suggesting that better knowledge of the cosmic ray 
source distribution and diffusion will have an important 
impa ct on the resul t s of 7 -ray studies. 

Ack ermann et al.] (|2012l ) carry out a thorough study 
of the impact of systematics on the global 7-ray analy- 
ses that include a cosmic-ray generation and propagation 
model, frequently used to infer Galactic Xqq gradients. 
They find that the Xco determination can be very sen- 
sitive to assumptions such as the cosmic ray source dis- 
tribution and the HI spin temperature, as well as the 
selection cuts in the templates (in particular, the dust 
template). Indeed, the value in the outer Galaxy in their 
analysis is extremely sensitive to the model cosmic ray 
source distribution. The magnitude of a Galactic Xco 
gradient turns out to also be very se nsitive to the un- 
derly ing assumptions (see Fig. 25 in lAckermann et all 
[2012). What appears a robust result is a uniformly low 
value of Xco near the Galactic center {R ^ — 1.5 kpc). 
This was already poin ted out by much earlier 7-ray stud- 
ies (|Blitz et al.l 119851) . The authors also conclude that 
including the dust information leads to an improvement 
in the agreements between models and 7-ray data, even 
on the Galactic plane, suggesting that "CO-faint" gas is 



an ubiquitous phenomenon. 

4.4. Synthesis: Value and Systematic Variations of 
Xco *^ the Milky Way 

There is an assuring degree of uniformity among the 
values of Xco obtained through the variety of method- 
ologies available in the Milky Way. Representative results 
from analyses using virial masses, CO isotopologues, dust 
extinction, dust emission, and diffuse j-ray radiation 
hover around a "typical" value for the disk of the Milky 
Way Xco ~ 2x10^0 cm-2(K km s^^)-! (Table[^. This 
fact, combined with the simple theoretical arguments 
outlined in ^ as to the physics behind the H2-to-CO 
conversion factor, as well as the results from elaborate 
numerical simulations discussed in ^3.2\ strongly suggests 
that we know the mass-to-light calibration for GMCs in 
the disk of the Milky Way to within ±0.3 dex certainly, 
and probably with an accuracy closer to ±0.1 dex (30%). 
This is an average number, valid over large scales. Indi- 
vidual GMCs will scatter around this value by a certain 
amount, and individual lines-of-sight will vary even more. 

There are, however, some systematic departures from 
this value in particular regimes, as suggested by the sim- 
ple theoretical arguments. As pointed out by several 
studies, the Galactic center region appears to have an 
Xco value 3 — 10 times lower than the di sk. In ad- 
dition to the a forementi oned 7-ray results (iBlitz et al.l 
119851 iStrong et al. 2004| lAckermann et al.ll2012t ). a low 
Xco in the center of the Milky Wa y was o btained by 
analysis of the dust emission ( Sodroski et a l.|fT9951. and 
the virial mass of its clouds (jOka et al.l ll998. 2001). We 
will see in ^ that this is not uncommonly observed in 
the centers of other galaxies. It is likely due to a com- 
bination of enhanced excitation (clouds are hotter) and 
the dynamical effects discussed in §2.31 

Departures may occur in some nearby, high-latitude 
clouds, for example the Cepheus/Polaris Flare region. A 
critical step to understand the magnitude or even ex- 
istence of such discrepancies will be to derive Xco for 
diffuse 7 ray observations using conventions matched 
to other fields. We have attempted some estimate 
of the correction in the Chameleon, R CrA, and the 
Cepheus/Polaris Flare. In Cepheus/Polaris even after 
trying to account for the difference in Xco convention, 
the 7-ray studies find a low Xco value suggesting that 
those local clouds are overluminous in CO by a factor 
^ 2, or maybe that their dust emissivity is different. Low 
values are not co nsistent with results from dust emis- 
sion techniques (iDame. Hartmann fc ThaddeusI 120011 : 
iPlanck Collaboration et al.ll2011a[) . On the other hand, 
they are known to be turb ulent to ver y small (~ 0.01 
pc) scale ()Mivi llc-Dcschcnc s et al.ll2010f ). a nd turbulent 
dissip ation may play a role in exciting CO (jlngalls et al.l 

[MK). 

Finally, observations of Xco in the outer galaxy are 
still sparse, although s ome of them suggest high values 
of Xco (for example, iHever. Carpenter &: Snelll 12001). 
This can be understood in terms of an increasing domi- 
nance of the "CO-faint" molecular gas associated prob- 
ably with decreasing metallicities in the outer disk (see, 
fJS]). In nearby, resolved molecular clouds there is com- 
pelling evidence that a regime of high Xco exists at low 
extinctions and column densities fe.g.. iGoldsmith et al.l 
i2008; .Pineda et al..,201Q) . where most of the carbon in 
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TABLE 1 

Representative Xco values in the Milky Way disk 



Method 


V /in20 


References 




cm (K km s } 




V irial 


O 1 


DOlonion ct al. ( iyo ( ) 




2.8 


_bcovillc ct al. (1987) 


Isotopologues 


1.8 


^Cjoidsmitn ct al. {zUUo} 


Extinction 


1.8 


l^rcrkmg, Langcr &: Wiison (lyoz) 






j^omoarQi. >\ives 6i ijaaa izuudi 




0.9 - 3.0 


Pineda. Caselli & Goodman (2008) 




2.1 


Pineda et al. (2010) 




1.7-2.3 


Paradis et al. ("2012) 


Dust Emission 


1.8 


Dame, Hartmann & Thaddeus f2001) 




2.5 


Planck Collaboration et al. (2011a) 


7-rays 


1.9 


Strong & Mattox ("1996) 




1.7 


Grenier. Oasandiian & Terrier (2005) 




0.9 - 1.9 * 


Abdo et al. (2010d) 




1.9 - 2.1 * 


Ackermann ct al. (2011, 2012d) 




0.7-1.0 * 


Ackcrmann ct al. (2012b ,£) 


* Note difference 


in Xco convention (i|4.3ll. 



the gas phase is not locked in CO molecules. 

It seems clear that CO becomes a poor tracer of H2 
at low column densities. Most studie s agree that this 
phen o menon occurs for Ay < 2 (iGoldsmith et al.l 
[2OOI IPlanck Collaboration etahl l2011aD . How 
much molecular gas exists in this regime is still a 
subject of study, but it seems to represent a sub- 
stantial, perhaps even dominant fraction of the H2 
near the s olar circle (iGrenier. Casandiian fc Terrieii 
2005t 'Wolfi re. HoUenbach fc McKe3 \2(M 

Planck Col laborati on et al. I l2011aD . It is, however, 
unlikely to be a large fraction of the total molecular 
mass of the Galaxy, which is dominated by clouds in 
the inner galaxy. Nonetheless, we will see that this 
"CO-faint" phase almost surely constitutes most of the 
molecular gas in low metallicity systems. 

Using these conclusions, what is the characteristic 
surface density of a CMC, Sgmc, in the Milky Way? 
The precise value, or even whether it is a well-defined 
quantity, is a matte r of study. Using the survey by 
iSolomon et al.l ()1987[) updated to the new distance scale 
and assuming Xco, 20 = 2, we find a distribution of sur- 
face densities with Sqmc ~ 150^70 Mq pc~^ (ilcr in- 
terval). Using instead the distribution of 331cr^/i? to 
evaluate surface density (see discussion after Eq. 
we obtain Egmc « 200^80° M© pc'^. IHever et al] 
(I2OO9I) use ^■'CO observations of these clouds, finding 
^GMC ~ 40 M0 pc~^ over the same cloud areas but 
concluding that it is likely an underestimate of at least 
a factor of 2 due to non-LTE and optical depth ef- 
fects (Sgmc ~ 144 M0 pc~^ using a more restrictive 
^•^CO cont our instead of the original ^^CO to define sur- 
face area, iRoman-Duval et all 120101) . Real cloud s have 
a range of surface densities. iHeiderman "eFall (l2010l ) 
analyze extinction-based measurements for 20 nearby 
clouds, calculating their distribution of Smoi- It is clear 
that many of these clouds do not reach the characteris- 
tic Egmc of clouds in the inner Galaxy; for the most 
part the clouds with low mass tend to have low sur- 
face densities. Nonetheless, several of the most mas- 



sive clouds do have Egmc ^ 100 M© pc ^. The mass- 
weighted (area-weighted) Egmc are 160, 150, 140, and 
110 (140, 140, 100, and 90) M© pc'^ for Serpens-Aquila, 
Serpens, Ophiuchus, and Perseus respectively according 
to their Table 2. By comparison, Egmc ~ 85 Mq pc~^ 
for a sample of nea rby galaxies, many of them dwarf 
()Bolatto et al.ll2008| ). Note that Egmc is likely to be a 
function of the envi r onment. In the Galactic Center, for 
example. lOka et aD ()1998L [200ll ) report a size-line width 
relation where the coefficient C is five times larger than 
that observed in the Milky Way disk suggesting a Egmc 
that is 25 times larger (see also iRosolowskv fc Blitd 
l2005[) . 

5. Xco IN NORMAL GALAXIES 

The distance of even the nearest galaxies renders CO 
the primary tracer of molecular gas outside the Milky 
Way, a situation that will improve but not reverse with 
ALMA. Other galaxies therefore represent the key appli- 
cation of ^co- They also offer a wider range of envi- 
ronments than the Milky Way and a simpler mapping 
between local ISM conditions, especially metallicity, and 
line of sight. As a result, for more than two decades ob- 
servations of the nearest galaxies have been used to test 
and extend calibrations of Xqq as a function of metallic- 
ity and other local ISM properties. Here we review the 
techniques available to derive Xco in "normal" galaxies, 
meaning star-forming dwarf, spiral, or elliptical galaxies. 
We discuss specific efforts to understand the behavior of 
^co as a function of metallicity in f|6]and the special case 
of overwhelmingly molecular starburst galaxies, such as 
the local luminous and ultraluminous infrared galaxies, 
in fJT] considers Xco in galaxies at high redshift. 

Only a subset of the techniques used to determine Xqo 
in the Milky Way can be applied to other galaxies. In 
each case, the limited sensitivity and resolution of mil- 
limeter, submillimctcr, and infrared facilities complicate 
the calculation. Direct estimates rely almost exclusively 
on the use of virial mass measurements f tj4.ip . dust emis- 
sion employed as an optically thin tracer of the total gas 
reservoir ( ii4.2p . or modeling of multiple CO lines. Ide- 
ally such modeling includes optically thin isotopologues. 
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but studies of high redshift systems must often make due 
with a few, or only one, ^^CO hue ratios (Sj8]). 

5.1. Extragalactic Virial Mass Estimates 

Since the late 1980s, millimeter telescopes have been 
able to resolve CO emission from nearby galaxies into 
discrete molecular clouds. From such observations, one 
can estimate the line width, size, and luminosity of these 
objects and proceed as in §4.11 (see the recent review by 
IFukui fc KawamurallMO ). 

Resolution and sensitivity have limited virial 
mass measurements to the nearest galaxies, those 
in the Local Group and its immediate environs. 
The Magellanic Clouds {d ~ 50 kpc) are close 
enough that sing le dish t elescopes resolve Individ- 



Israel et all 20031 iMizuno et al. 


120011: IFukui et 3l 


20081 IHuEhes et al.l 20101 


Wong et al.l 120111). 



Millimeter-wave interferometers resolve populations 
of GMCs in other Local Group galaxies — M 31, 
M 33, NGC 6822, and IC 10 (d ^ 1 Mpc e.g., 
Vogel. Boulanger fc Ba lll \l98% IWilson fc Scoville TOOOI: 
Wilson 1995; Rosolow skv et al.lT2 003: J^rov et al. 200l 
Blitz et al.ll2007t IFukui fc Kawamurall2010[ ). 



Early measurements beyond the Local G roup focused 
on ve r y nearby dwa r f gala. xi es (d ~ 3 Mpc. iTavlor et al.l 
ITQqI iWalter et al] [200ll [200l IBolatto et al.l 120081 ) 
or considered structures more massive than GMCs 
in more distant gala xies (Giant Molecular Associa- 
tions, or GMAs, e.g., iVogel. Kulkarni fc Scovilld 119881: 
iRand fc Kulkarnilll990[ ). The nearest spiral galaxies tend 
to be more distant (d ~ 6 Mpc), so that ~ 1" reso- 
lution is required to resolve massive GMCs. Only re- 
cently mm-wave interferometers have begun to achieve 
this resolution with the requisite sensitivity, by invest- 
ing large amounts of time into dedicated observations 
of b right regions of the neares t mass i ve spiral galax- 
ies (iDonovaii Mever et all 120121 120131: iRebolledo eFal] 
12012 ). This has allowed the first cloud-scale virial mass 
measurements of populations of clouds in spiral galaxies 
beyond the Local Group. 

5.1.1. Xco From Extragalactic Virial Mass Analyses 

Broadly, virial mass measurements across a wide range 
of environments yield Xco,20 ~ 1 — 4, consistent with 
Milky Way results and very similar to Xqo derived from 
dust -based techniq u es app lied to high mass nearby galax- 
ies. IBolatto et aP (pOOl find Xco,20 « ^-^tH in 12 
nearby systems. In the highest resolutio n studies of 
the La r ge Magellanic Cloud (L MC) to date , Israel et all 
(I2OOI1. IHu ehes et al.l ((20Tol) . IWong et al.l ^01^, and 
I Pineda et al. (201(|)all find Xco,20 ~ 4. Considering 
NGC 6946, on e of the near e st mo lecule-rich spiral galax- 
ies, IPonovan Mever et al.l ()2012f ) report Xqo 20 ~ 1-2 
and IRebolledo et al.l (|20T2|) find Xcoao ~ 2. 'in M 33, 
iRosolowskv et al.l ( 20031) show -^00. 20 ^ 2 independent 
of radius and metallicity. In M 31, iRosolowskvl ()2007f ) 
calculate Xco.2a ~ 4 assuming virial equilibrium. 

Figure [6] shows the relationship between virial mass 
and CO luminosity for a subset of these measurements. 
A good correlation extends across three orders magni- 
tude in luminosity and roughly a dozen systems. As in 
the Milky way, there may be evidence for a slightly sub- 
linear slope (see ^2. II and ^^4. II and Figure [2]), though this 
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Fig. 6. — Relation between virial mass (ly-axis) and CO lumi- 
nosity (x-axis) for extragalactic GMCs. We show virial masses 
measured from selected high spatial resolution CO o bserv ations of 
nearby galaxies: the compilation of I Bolatto et al] l|2008l . includ- 
ing M31, M33, an d nine dwarf galaxies), high reso lutio n stud- 
ies of the LMC bv IPineda et al.l ( 2009) and Wong clTaLl l|20Tll V 
and high resolution studies of the nearby spiral N GC 6946 by 
IPonovan Mever et al.l | |20T^ and iRebolledo et al. 1201^'), as we ll 
as NGC 4826 and NGC 4736 bv IPonovan Mever et al.l (120131 ). 
Pashed lines show fixed Xqo^ with the typical Milky Way value 
Xco,20 = 2 and ±30% indicated by the gray region. Virial mass 
correlates with luminosity, albeit with large scatter, across more 
than three orders of magnitude in extragalactic systems. The me- 
dian across all displayed data is Jfco,20 = 2.8 with 0.4 dex scatter 
and the best fit relation has a power law index 0.90 it 0.05, reflect- 
ing that most Local Group clouds show Xqq 20 ~ 4 while both 
studies of the bright spiral NGC 6946 find a lower Xco,20 ~ 1 — 2. 

may also reflect methodological or environmental differ- 
ences among the galaxies studied. Across the ensemble 
of points we plot the median Xco.20 = 2.9 with 0.4 dex 
scatter (slightly larger than a factor of 2). 

Surprisingly at first, these "Galactic" Xco values 
obtained from virial masses extend to low metallic- 
ity, irregular galaxies. In the low metallicity Small 
Magellanic Cloud (SMC, wit h 12 -HogfO/H ] » 8.0 , 
Duf our. Shields fc TalbotJ fT982l ). IBolatto et all (I200I 
and ilsrael et al.l ( 20031 ) find Xcn.2ri » 4, c onsist ent with 
the Milky Way value. IRosolowskv et al.l (I2003D find no 
depen dence of Xco on metallicity in M33. iLerov et ahl 
(|2006f ) find an approximately Galactic X cp in the Local 
Grou p dwarf IC 10 (1 2 + log[0/H] w 8.17. lLeaueux et all 
[Ml). IWilsonI (p9l f ound Xr.c, 20 < 6.6 in NGC 6822 
(12 + log[0/H] R:! 8.20, ILeaueux et all HMp). Most of 
the ga laxies studied by the aforementioned lBolatto et all 
((2OOI are dwarf irregulars with subsolar metallicity. In 
each case we highlight the results for the highest reso- 
lution study of the galaxy in question. We discuss the 
effects of varying spatial resolution in the next section. 
Overall, these studies show that high spatial resolution 
virial mass measurements suggest roughly Galactic Xqo 
irrespective of metallicity. 
Virial masses have also been measured for 

GMAs or superclouds 



very 



masses 
large structures, 
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(|Vogel. Kulkarni fc Scovill^ ITOSSl: IWilson et alJ [200l . 
although it is unclear t he degree to which they are 
virialized or even bound. IWilson et al.l ()2003l ) consider 
^ 500 pc-scale structures in the Antennae Galaxies 
and find the Xng to be a pproximately Galactic (see 
also lUeda et aljl2012l : iWei. Ke to & Ho 2012) . but this 
conflicts with results from spectral line modeling (see 
^ . Phrasing their resu l ts larg ely in terms of bounded- 
ness, iRand fc Kulkarnil (|1990f) and lAdler et alJ ()1992f ) 
find Xco,20 ~ 3 and Xco,20 ~ 1-2 to be needed for 
yirialized GMAs in M51. Studies of other ga laxies, for 
example M 83 ()Rand. Lord fc Higdonl Il999[ ). suggest 
that not all GMAs are grayitationally bound and that 
the mass spectrum of GMAs yaries systematically from 
galaxy to galaxy o r betw e en arm and interarm regions 
Mand fc Kulkarni [199(11 : IRand. Lord fc Higdonl [l999l : 
IWilson et al.ll2003f ). complicating the interpretation of 
these large-scale measurements. 

5.1.2. Caveats on Virial Mass-Based Xco Estimates 



Wong et al 



Fukui et al 



As we will see in i )5.21 the approximately constant 
Galactic Xqq implied by virial masses on small scales 
at low metallicity appears to contradict the finding from 
dust-based measurements and other scaling arguments 
(for example iBlanc et al.l 120131) , which consistently in- 
dicate that Xco increases with decreasing metallicity. 
This discrepancy most likely arises because virial mass 
measurements sample the gas that is bright in CO, while 
dust-based measurements include all H2 along the line-of- 
sight. As discussed below (S|6|), decreasing dust shielding 
at low metallicities causes CO to be preferentially pho- 
todissociated relative to H2 , creating a massive reservoir 
of H2 in which C+ and C rather than CO represent the 
dominant forms of gas-phase carbon. Because this reser- 
voir is external to the CO emitting surface it will not be 
reflected in the CO size or its line width (unless the sur- 
face pressure term is important). Thus, we expect that 
high resolution virial mass measurements preferentially 
probe Xco in the CO-bright region. 

The range of values discussed above, Arco,20 ^1 — 4, 
is significant. Does it indicate real variations in Xco? 
Measuring cloud properties involves several methodolog- 
ical choices f ij4.ip . Different methods applied to the 
same Milky Way data shift results by ~ 30%, and bi- 
ases of ^ 40% are commo n in extragalactic data (e.g., 
iRosolowsky fc Lerovl 120061 ) . As we discuss below, the 
A"co obtained by virial mass analysis seems to depend 
on the physical resolution of the observations. Ideally, 
results will be compared to "control" measurements that 
have been extracted and analyzed in an identical way, 
ideally at matched spectral and spatial resolution and 
sensitivity (e.g., Pineda et al. 2009). Many studies now 
empl oy the CPROPS algorithm (Rosolows kv fc Lerovl 
l2006f ). which is designed to account for sensitivity and 
resolution biases in a systematic way, allowing ready 
cross-comparison among data sets. This is not a perfect 
substitute, however, for matched analyses. Our assess- 
ment is that in lieu of such careful comparison, differences 
of ±50% in Xco should still be viewed as qualitatively 
similar. 

We raised the issue of spatial scale in the discus- 
sion of systematics and implicitly in the discussion 
of GMAs. Virial mass-based Xco exhibit a complex 
dependence on the spatial scale of the observations. 



For dwarf galaxies it seems that even using similar 
methodology and accounting for resolution biases, stud- 
ies with finer spatial resolution systematically return 
low er Xco than coarser resolution studies. For exam- 
ple, 'Hugh es et all (IMol) find Xco m 4 in the LMC 
(as do .Israel et al.l 120031 : IPTneda et al.ll2009l : 

m compared to v^rco.ao ~ 7 found by 
using similar methodology but wi th ^ 3 times 
coarse r linear resolution. In the SMC, iMizuno et al.l 
(|2001h find Xco, 20 ~ 14, while higher resolution stud- 
ies find Xco as low as Xcq^jq^ 2 — 4 for the small- 
est r esolved objects ([Israel et al] [20031 : iBolatto et al.l 
20031). Contra sting the interferometer measurements of 
Wilson! (I1994D and the coarser single-dish observations 
of iGratier et al.l (|2010D reveals a similar discrepancy in 
NGC 6822. In one of the first stud i es to consider this 
effect, 'Rubio. Lequeux fc Boulanged (jl993) explicitly fit 
a dependence for Xco in the SMC as a function of spa- 
tial scale, finding Xco <x R°' ^ (see also the m ultiscale 
analyses in IBolatto et a"ni2003l : iLerov et al.|[2009ll . 

This scale dependence may reflect one of several sce- 
narios. First, low resolution observations can associate 
physically distinct clouds that are not bound, causing a 
virial mass analysis to overpredict Xco- Alternatively, 
the ensemble of clouds conflated by a coarser beam may 
indeed be bound. In the case of a heavily molecular in- 
terstellar medium like the Antennae Galaxies or the arm 
regions of M 51, most of the material in the larger bound 
structure may be molecular and a virial mass measure- 
ment may yield a meaningful, nearly Galactic conversion 
factor for objects much bigger than standard GMCs. For 
a low metallicity irregular galaxy, the best case for low 
resolution virial mass measurements is that large com- 
plexes are virialized and that the low Ay gas between 
the bright clouds is H2 associated with [CII]. We caution, 
however, that this is only one of many possible scenarios. 

We suggest that the sensitivity of virial mass mea- 
surements to extended CO-free envelopes of H2 is am- 
biguous at best. In that sense, the uniformity in Xco 
derived from virial masses probably reflects fairly uni- 
form conditions in C O-bright regions of molecular clouds 
(jBolatto et al.l[2008h . The spatial scale at which virial- 
ized structures emerge in galaxies is unclear. Given these 
ambiguities and the assumptions involved in the calcula- 
tion of virial masses, we emphasize the need for careful 
comparison to matched data to interpret virial results. 



5.2. Extragalactic Dust-Based Estimates of Xco 

Dust is expected and observed to be well-mixed with 
gas, and dust emission remains optically thin over most 
regions of normal galaxies. Following the approach out- 
lined in §4.2| dust emission offers a tool to estimate Xco ■ 
Modeling infrared or millimeter emission yields an esti- 
mate of the dust optical depth, r^. Using Eq. [221 Xco 
can be estimated from Td, HI, CO, and the dust emissiv- 
ity per H atom (Eq. [22]) . Conventions in the literature 
vary, with sometimes combined with a dust mass ab- 
sorption coefficient and used as a dust mass, and <5dgr 
alternately cast as the emissivity per H atom or a dust- 
to-gas mass ratio. Regardless of convention, the critical 
elements are a linear tracer of the dust surface density 
and a calibration of the relation between this tracer and 
gas column density. 
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Fig. 7. — Dust-based aco de termination across the disks of 22 
galaxies fSandstrom et aTll2012l ). The top panel shows the results 
for the individual regions (gray points), averages in radial bins for 
each galaxy (color lines), and the average trend for all data as a 
function of radius (black circles with error bars showing the scatter 
in that bin). The black solid line in the top panel corresponds to 
^CO,20 = 2, and the dotted line shows the average, weighting 
all solution pixels equally (Xco,20 = 1.2). Weighting by galaxy 
or CO intensity rather than line of sight, the mean is Xqq 20 ~ 
1.4 — 1.8 for low inclination galaxies. The bottom panel shows 
each galaxy normalized by its mean Xqq to highlight the fact that 
many galaxies exhibit low ^co ™ their centers. 



iThronson efal] (|1988D and lThronsoU (|1988bl ) first sug- 
gest and apply varia tions on this technique to nearby 
galaxies. Ilsraeil ()1997D uses data from IRAS to carry out 
the first comprehensive dust-based extragalactic Xqo 
study. He considers individual regions in eight (mostly 
irregular) galaxies and derives from a combination of 
60 and 100/im continuum data. [Israeli estimates fccR: 
the dust emissivity per H atom, from comparison of 
and HI in regions within the galaxy of interest but cho- 
sen to lie well away from areas of active star formation, 
thus assumed to be atomic-dominated. This internally 
derived estimate of fccR represents a key strength of the 
approach: Sbgr is derived self-consistently from compar- 
ison of HI and r^, and does not rely on assuming a dust- 
to-gas ratio. This leads to the cancellation of many sys - 
tematic errors (e.g., see llsraeilll997bl : iLerov et al.ll2007D . 
leaving only variations in (5dgr and within the target 
to affect the determination of Xco in Eq. [231 

5.2.1. Dust-Based Estimates m Normal Disk Galaxies 

Spitzer and Herschel allowed the extension of the dust 
approach to more mass ive, mor e dista nt, and more "nor- 
mal" galaxies. iDraine et al.l (|2Q07D compare galaxy- 
integrated infrared spectral energy distribution (SED) 
modeling to CO and HI luminosities for a large sample. 
They argue that Xco, 20 ~ 4 (over entire galaxies) yields 
the most sensible gas-to-dust ratio results in their sam- 
ple. ILerov et al.l (|2011) perform a self-consistent treat- 
ment of the Local Group galaxies M 31, M 33, LMC, 
NGC 6822, and SMC. They find Xco,20 « 1 - 4.5 for re- 
gions of M 31, M 33, an d the LMC (the high er metallicity 
galaxies in the sample) . iSmith et al.l (|2012D use Herschel 
observations to solve for Xco finding Xco, 20 ~ 2 in 



M 31. 

Other dust tracers such as mm- wave continuum emis- 
sion and visual extinction have been used to arrive at 
Xco estimates, although usually by assuming or scal- 
ing a Galactic calibrat ion with the associated system- 
atic uncertainties. .Guelin et al.l f|1993) use 1.2 mm con- 
tinuum observations of the edge-on spiral NGC 8 91 to 
estimate Xc o .2n 1- Both iNakai fc Kunol ()1995[) and 
iGuelin et all (|1995D study Xco in M 51, the first using 
extinction estimates from Ha/H/3, the second millimeter- 
wave continuum emission. Nakai fc Kun o ( 199S ) arrive 
at Xco, 20 ~ 0.9 ±0.1 with a factor of ^ 2 variatio n with 
galactocentric radius, wh i le IGuelin et all ()1995[ ) finds 
Xco,2o ~ 0.6. iZhu et al.l ([20091) employs 850 /xm data 
to check the Xco derived from spectral line modeling in 
NGC 3310 and NGC 157, finding good agreement with 
Xco. 20 ~ 1-3 — 3.0 in the disks of their targets and a 
much lower Xco in the ce nter of NGC 157. 

iSandstrom et al.l (|2012D carry out the most compre- 
hensive extragalactic study of Xco to date. They com- 
bine high-quality CO J = 2 — > 1 maps with Herschel 
and Spitzer dust continuu m, and h igh-resolution HI data. 
Building on the method of ILerov e t al. (2011), they break 
apart galaxies into regions several kpc^ in size and within 
each region they simultaneously solve for (5dgr and Xco • 
This yields resolved, self-consistent Xco measurements 
across 22 galaxy disks. Their methodology requires good 
S/N CO detections and so restricts robust Xco measure- 
ments to reasonably CO-bright parts of galaxies, typi- 
cally half of the optical disks. There, the authors find 
Xco, 20 ~ 1.4-1.813 with a la scatter among individual 
solutions of 0.4 dex. 

The top pane l of Fig. 13 shows Xco derived by 
ISandstrom et al.l as a function of galactocentric radius 
for their whole sample. The bottom panel shows Xco 
for each region normalized to the average value for the 
galaxy, highlighting the internal radial structure. This 
study finds clear central Xco depressions for a number 
galaxies. The value of Xco in galaxy centers relative to 
their disks spans a range from no change to a factor of 
five below the galaxy average. The authors note that 
these low central Xco values tend to coincide with high 
stellar masses and enhanced CO brightness, suggesting 
that the effects discussed for starbursts in ^Ilniay also be 
at play in the central parts of many galaxies, including 
the Milky Way 

Thus in the disks of normal star-forming galaxies, a 
dust-based approach yields Xco, 20 ~ 1 — 4 on kpc 
scales. For comparison, in the Milky Way extinction 
yields Xco 20 ~ 1-7 — 2.3 and dust emission yields 
Xco,20 ~ 1-8 - 2.5 (S21 Table H]). Thus, broadly dust 
analyses strongly support a Milky Way conversion factor 
in the disks of normal, massive disk galaxies but method- 
ological differences and, presumably, real changes in Xco 
with environment produce a factor of 2 spread among 
studies and galaxies. 

5.2.2. Dust-Based Estimates in Dwarf Irregular Galaxies 



^ The weighting used to derive the average Xqq affects 
the mean value. Weighting all high quality solutions equally, 
ISandstrom et aT] find mean Xco, 20 ~ 1-4, median Xco, 20 ~ 1-2- 
Weighting instead by CO intensity, they find Xco, 20 ~ 1-8. 
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Dust-based determinations of Xqq in low metal- 
licity dwarf irregular galaxies consistently yield val- 
ues much higher than Galacti c, an d also generally 
higher than virial estimates. Ilsraell ()1997f ) finds in 
the SMC, the lowest meta llicity target, a notably high 
-^^00,20 ~ 120 ± 20. .Boselli. Lequeux fc Gava zzi (2002) 
carry out a similar calculation and also find system- 
atically higher Xqo at low metallicity. Subsequent 
studies leveraging the sensitivity, resolution, and wave- 
length coverage of Spitzer and Herschel qualitatively sup- 
port t his co nclusion. [L croy ct al. (2007), Lcro v et aTl 
(l2009l). and IBolatto et all (1201 ID follow risrael (|1997[ ) 
and Istanimirovic et al.l ( 2000l ) in analyzing the SMC. 
They find Xco,20 « 40 — 120 using measurements 
that isolate in dividual clouds and estimate (5dgr locally. 
IGratier et all (2010) considers NGC 6822 using a sim- 
ilar technique, also finding high Xco,20 ~ 40 on large 
scales. Attempting to minimize systematics a nd m aking 
very conservative assumptions, iLerov et"al] ()2011f ) still 
find Xco,20 « 10 and 40 for NGC 6822 and the SMC, 
significantly higher than they find for the higher metal- 
hcity M 31, M 33, and the LMC. 

In the Magellanic Clouds, Herschel and millimeter- 
wave bolo meter cameras have mapped individual molecu- 
lar clouds. iRoman-Duval et al.l (|2010f ) find an increase in 
Xco in the poorly s hi elded oute r region s of L MC clouds. 
iRubio et al.f (|200l . iBot et al.l (|200l . and iBot et al.l 
(|2010t ) find une xpectedly bright m illimeter-wave emission 
in SMC clouds. IBot et aTl ([2010') show that the emission 
suggests a factor of four higher masses than returned 
by a virial analysis, even on very small scales and even 
assuming a very conservative rescaled Milky Way dust 
mass emissivity coefficient. 

Therefore, dust-based Xqq estimates indicate high con- 
version factors in low metallicity, irregular systems. This 
agrees with some other methods of estimating Xqq but 
not with virial mass results, particularly on small scales, 
as already discussed (t jS.ip . We discuss the discrepancy 
further below (S}6]). The simplest explanation is that 
virial masses based on CO emission do not sample the 
full potential well of the cloud in low metallicity systems, 
where CO is selectively photodissociated relative to H2 
at low extinctions. 

5.2.3. Caveats on Dust-Based Xco Estimates 

Because dust traces the total gas in the system with 
little bias (at least compared to molecules), this tech- 
nique represents potentially the most direct way to esti- 
mate Xco in other galaxies. We highlight two important 
caveats: the possibility of other "invisible" gas compo- 
nents and variations in the emissivity. 

As described, the method will trace all gas not ac- 
counted for by HI 21 cm observations and assign it to 
H2. The good qualitative and quantitative association 
between H2 derived in this manner and CO emission sug- 
gests that H2 does represent the dominan t component 
(e.g., iDame. Hartmann fc ThaddeusI 120011 : ILerov et"al] 
|2009| ). but opaque HI and ionized gas may still host 
dust that masquerades as "H2" using Eq. [23l Both 
components could potentially represent significant mass 
reservoirs. Their impact on dust-based A^co estimates 
depends on the methodology. Any smooth component 
evenly mixed with the rest of the ISM will "divide out" 
of a self-consistent analysis. Moreover, the warm ion- 



ized medium exhibits a large scale height compared to 
the cooler gas discussed here and may be sub ject to 
dust d estruction without replenish i nent (e .g., see f Praiiii 
120091 ) . iPlanck Collaboration etHI (|2011bl ) find that dust 
associated with their ionized gas template represents a 
minor source of emission. On the other hand, absorp- 
tion line experiments from our galaxy and others sug- 
,est that ^ 20% H I may be missed due to opacity effects 
Heiles fc Trolandl 12003) and observations of M 31 sug- 
gest that it may lie in a filamentary, dense component 
(jBraun et a l. 2009), making it perhaps a more likely con- 
taminant in Xco determinations. Nonetheless, the close 
agreement between Xco derived from dust and other de- 
terminations in our Galaxy strongly suggests that con- 
tamination from dust mixed with other "invisible" com- 
ponents, particularly the ionized gas, represents a minor 
correction. 

Concerning the second caveat, the methodology relies 
on either blindly assuming or self-consistently determin- 
ing the dust emissivity per hydrogen, Jdgr- This is a 
combination of the dust-to-gas mass ratio, and the mass 
absorption coefficient of dust at IR wavelengths. In blind 
determinations, a linear scaling of Eq. [22] with metallic- 
ity is usually assumed. In self-consistent Xco determi- 
nations, 5dgr is measured locally through the ratio of 
the dust tracer to HI somewhere within the galaxy (for 
example, in an Hl-dominated region), and ideally close 
to or inside the region where Aco is estimated. This 
makes self-consistent determinations very robust. Even 
the most robust analyses, however, still assume (5dgr 
to be constant over some region of a galaxy and along 
a line of sight, and between the atomic and molecular 
phases. The (5dgRj however, may change across galaxies 
and between phases due to metallicity gradients, varying 
balance of dust creation and destruction, or changes in 
dust grain proper ties. 

As discussed in ^ ^4.21 recent Planck results suggest only 
mild localized variations in the emissivity per H nucleon 
across the Milky Way. Observations of other g alaxies 
do show that (5nGR. depends on m etallicity (Drain eet al.l 
I2007t iMufioz-Mateos et"alll2009t ). Large scale ISM den- 
sity may also be important; (5dgr appears depressed 
in the low-density SMC Wing and the outer env elopes 
of dwarf galaxies (|Bot et all [20041; ILerov et al.l 120071: 
iDraine et al.ll2007|) . The balance of the dust production 
and destruction mechanisms is complex: basic account- 
ing i mplies that most dust ma ss buildup occurs in the 
ISM (jDwekl 119981 : lDraindl200"9l and references therein), 
suggesting an increase in (5dgr in denser regions of the 
ISM where this accumulation must take place. Inde- 
pendent of dust mass buildup, if grains become more 
efficient emitters for their mass at moderate densities 
(through the formation of fluffy aggregates, for exam- 
ple) then the effective Sdcr will be higher at high densi- 
ties. Though not strongly favored by the overall Planck 
results, which suggest that emissivity variations are lo- 
calized, a number of authors point to evidence for an 
environmental dependency of the dust emissivity {e.g., 
Cambrcsy, Jarrctt & Bcichman 2005; Bot et al. 2003). 

In our opinion these caveats concerning emissivity will 
lead a dust-based approach to preferentially overpredict 
the amount of H2 present and consequently Aco 1 as dust 
associated with the molecular ISM will have a higher 
dust-to-gas ratio or be better at emitting in the far- 
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infrared. Nonetheless, we still see the dust-based ap- 
proach as the most reliable way of producing extragalac- 
tic Xco estimates. Placing stronger quantitative con- 
straints on these systematics requires further work. The 
Planck results noted above, and the close agreement in 
the Milky Way of dust-based Xco with other techniques, 
suggests that emissivity variations are not a major con- 
cern. 

5.3. Extragalactic Spectral Line Modeling 

Observations of multiple CO lines or a combination of 
CO and other chemical species allow one to constrain 
the physical conditions that give rise to CO emission. 
When these observations include optically thinner tracers 
like ^'^CO, these constraints can be particularly powerful. 
Due to sensitivity considerations, most multi-line data 
sets have been assembled for bright regions like galaxy 
centers (or starbursts, SJT]). From these data we have 
constraints on Xqq in bright regions with very different 
systematics than either virial or dust-based techniques. 

Generally, these results indicate that Xco lower than 
Xco, 20 = 2 is common but not ubiquitous in the 
bright, central regions of galaxies. T his is in qualita- 
tive a greement with the independent iSandstrom et afl 
(I2012D dust-based results a l ready mentioned. For exam- 
ple. [Israel. Tilanus fc Baaj ()2006[) use large velocity gra- 
dient (LVG) modeling of multiple ^^CO, "CO, and [Cl] 
lines in the central region of M 51. They find Xco, 20 ~ 
0.25 — 0.75 in the cen tral regions, in good agreement 
with ear her w o rk by iGarc ia-Burillo. Comb es fc GerinI 
(|1993| ). [Israeli (|2009t ) and .Israel (2009b3 extend this 
work to the centers of ten bright and star burst galax- 
ies, finding Xco 20 ""^ -1 — 0.3 (see also 'Israe l et al.l 
11)03; Israel &J3aa|[200i|), again significantly lower than 
Galactic. B ased on maps of the ^^C O/^^CO ratio across 
NGC 3627, IWatanabe etldl (I20TI argue for a similar 
central depression in Xco, with dynamical effects asso- 
ciated with a galactic bar leading to broader line widths 
and optically thinner ^^CO emission. Lower conver- 
sion factors are not a un iversal result, however. Utiliz- 
ing high resolution data, iSchinnerer et al.l (|2010l ) found 
Xco much closer to Galactic in the arms of M 51, 
Xco, 20 ~ 1-3 — 2.0. They attribute the difference with 
previous studies to the implicit emphasis on GMCs in 
their high resolution data, speculating that it removes a 
diffuse CO-bright component that would drive Xco to 
lower values. 

In addition to degeneracy inherent in the modeling, 
these line ratio-based techniques suffer from the funda- 
mental bias of the virial mass technique. They are only 
sensitive to regions where CO is bright and so may miss 
any component of "CO-faint" II2. A handful of observa- 
tions of low metallicity dwarf galaxies have attempted to 
address this directly by combining CO and [CII] observa- 
tions, with [CII] employed as a trace r of the "CO-faint" 
mole cu lar regime iM aloncv & Black 119881 : IStacev et al.l 
Il991| ). iPoglitsch et al.l ([1995) find that the [ClI]-to-CO 
ratio in the 30 Doradus region of the LMC is ^ 60, 000, 
roughly an order of magnitude higher than is observed 
in Milky Wa y star formiiig regi o ns or other s t ar-for ming 
spirals (e.g..lStaceY et al.lll99lD . Ilsrael et al.l () 19961 ) and 



Ilsrael fc MalonevI (j201in " also find high values elsewhere 
in the LMC and SMC, with large scatter in the ra- 



tio. iMadden et al.l ()1997D find that regions in the Local 
Group dwarf IC 10 also exhibit very high [CII] 158^m 
emission compared to CO, 2-10 times that found in the 
Milky Way. They argue that [CII] emission cannot be 
readily explained by the available ionized or HI gas, and 
is thus an ind i cator of H2 where CO is photodissociated. 
iHunter et al.l ()200lD present ISO measurements for sev- 
eral more dwarf irregular galaxies, showing that at least 
three of these s ystems also exh i bit ve ry high [ClI]-to-CO 
emission ratios. 'Cormi er et al.l ([2010') present a Herschel 
map of the 158/im [CII] line in NGC 4214, again show- 
ing very high [ClI]-to-CO ratios (- 20, 000-70, 000) that 
cannot be explained by emission from [CII] associated 
with ionized gas. Inferring physical conditions from the 
[CII] line requires modeling. When such calculations are 
carried out generally imply a massive layer on the pho- 
todissociation region in which the dominant form of hy- 
drogen is II2 , while carbon remains ionized as C"*" (e.g., 
IMadden et alllT997l: IPak et al.llT99l . 

Thus, detailed spectral line studies of CO-bright sources 
often, but not always, suggest lower Xco the bright 
central regions of galaxies. Meanwhile, at low metallici- 
ties, [CII] observations suggest an important reservoir of 
H2 not traced by CO. 

5.4. Synthesis: Xco Normal Galaxies 

Taken together, the picture offered by extragalactic 
Xco determinations in normal galaxy disks resembles 
that in the Milky Way writ large ( §4.4p . Virial masses, 
dust-based estimates, and spectral line modeling all sug- 
gest Xco, 20 ~ 1 — 4 in the disks of normal spiral galax- 
ies. Systematics clearly still affect each determination 
at the 50% level, with physical effects likely adding to 
produce the factor-of-two dispersion. Given this, apply- 
ing a "Milky Way" Xco = 2 x lO^o cm-2(K km s^i)"! 
with an uncertainty of « 0.3 dex appears a good first- 
brush approach for normal star-forming galaxies. This 
applies to galaxies where the CO emission is dominated 
by self-gravitating clouds or cloud complexes with masses 
dominated by H2. 

Several very strong lines of evidence, as well as simple 
arguments, show that Xco increases sharply in systems 
with metallicities be low 12 -|- log[0/II] » 8.4 (approxi- 
mately one- half solar. lAsplund et al.ll20Q9l ). High spatial 
resolution virial estimates find approximately "Galactic" 
Xco in low metallicity clouds, but they are not sensi- 
tive to extended "CO-faint" H2 envelopes, only to CO- 
bright regions. We expand on the dependence of Xco 
on metallicity in ^ We caution about the usefulness 
of virial estimates on large spatial scales, particularly in 
low metallicity galaxies. But even in normal galaxies, 
concerns exist about whether molecular cloud complexes 
and associations are bound and dominated by H2 on large 
scales. 

Finally, in the central parts of galaxies spectral line 
modeling suggests that Xco is often, but not always, de- 
pressed in a manner similar to that seen in the Galactic 
center and in molecule-rich starbursts. Dust observations 
are consistent with this picture, revealing central depres- 
sions in Xco in some galaxies. Broader line widths, in- 
creased excitation, and the emergence of a diffuse molec- 
ular medium likely contribute to more CO emission per 
unit mass, decreasing Xco- We expand on this topic 
when we consider starbursts in ^ 
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6. Xco AT LOW METALLICITIES 

CO emission is frequently very faint or nonexis- 
tent in gas-rich, actively star-forming low metallicity 
galaxies. This is true in absolute luminosity terms, 
with CO emission almost completely absent in low 
mass dwarf irregular galaxies below metallicity 12 -t- 
log 0/H » 8.0 (Elmcgrccn, Morris & Elmcgrcc ni 119801: 
Tacconi fc Younel Il987t ITavlor. Kobulnickv fc SkillmanI 
1998( 1. Even when it is detected, as in the case of 
the very nearby S mall Magellanic Cloud, CO is faint 
(jlsrael et al.lll986[) . CO is also faint in a normalized 
sense. The CO-to-FIR ratios in blue, vigorously star- 
forming dwarf galaxies are clearly low er than in sp i- 
ral galaxies (Tacconi & Young 1987; You ng et al.lll996l ). 
and the same ap plies to other star formation rate trac- 
ers such as Ha (jYoung et al.l 119961 ). The CO per unit 
starli ght is also depressed in very low mass dwarfs 
(e.g.. lYoung fc Scovilld |1991[ ). and the ratio of CO to 
HI is markedly decreas ing for later morphological type 
(|Young fc Knczck 19891 

Very importantly for CO emission, the low metal abun- 
dance in these systems implies lower C and O abun- 
danc e s and low dust-to-gas ra tios (e.g., IDraine et al.l 
I2007t iMuiioz-Mateos et al.ll2009l ). Both of these factors 
will exert a strong effect on the relative distributions of 
H2 and CO. Dust serves as the site of H2 formation and 
also provides much of the far-ultraviolet shielding nec- 
essary to prevent molecules that are not strongly self- 
shielding, such as CO, from photodissociating. Note that 
galaxy mass is strongly correlated not only with metallic- 
ity, but also with a number of other parameters that may 
affect the equilibrium abundance of H2, such as spiral 
density waves and interstellar pressure, further compli- 
cating the picture of the relation between CO luminosity 
and molecular mass. The critical question is whether low 
CO luminosities indicate a true deficit of H2, or merely 
the suppression of CO emission in low metallicity envi- 
ronments. 

Theoretical considerations, described below, lead us to 
expect a change in Xqq at low metallicities. Several 
strong lines of observational evidence, described in SJSI 
also suggest that Xco increases sharply in irregular, low- 
metallicity systems. Dust-based Xqo estimates, which 
are sensitive to the total gas present, return high values in 
such targets. The [ClI]-to-CO ratio seems dramatically 
higher in dwarf irregular galaxies than in spiral galaxies, 
supporting the hypothesis of a reservoir of H2 at low Ay 
where most carbon is locked in [CII] rather than CO. The 
SFR-to-CO ratio also increases, offering circumstantial 
evidence for star-forming gas not traced by CO. 

6.1. Theoretical Expectations for Low Metallicity Gas 

Naively, one would expect that a lower C and O would 
imply consequently fainter CO emission. Because CO 
is optically thick, however, its luminosity is determined 
by the emitting area and its brightness temperature and 
velocity spread. The dependence of Xqq on metallic- 
ity hinges fundamentally on how the relative sizes of the 
[CII] and CO-emitting regions change with lower heavy 
element and dust abundan ce (Fig. |Hl see the discu ssion 
of cloud structure in ^T^ . iMalonev fc Blackl ()1988D car- 
ried out one of the first thorough analyses of the prob- 
lem. Here we note the physical drivers for the location 



of bright CO emission and then discuss the implications 
for low metallicity. 

The relative distributions of CO and H2 will be a de- 
tailed function of the balance of formation and destruc- 
tion of C O, processes dis cussed for diffuse and dense ma- 
terial by Ivan Dishoeck fc Black (1986, 1988) (see also 
iWolfire. Hollenbach fc M cKcc 2010, appendices B and 
C). CO formation proceeds mainly through the produc- 
tion of OH via ion-neutral reactions initiated by cosmic- 
ray ionization. Once OH is formed the chemistry pro- 
ceeds via ion-ncutral reactions to form HCO+ followed 
by the dissociative recombination of HCO+ to form CO. 
The formation rates of CO thus depend on the abun- 
dance of C times the abundance of O. Because the rate 
of formation depends on the abundance of OH, the rate of 
destruction of OH by far-ultraviolet photons represents 
an additional limiting factor in setting the CO abundance 
(e.g.. .Wolfire. Hollenbach fc McKee..20ia Appendix C). 




Fig. 8. — Effect of metallicity on CO and If 2 in a spherical 
clump immersed in a uniform radiation field. Blue shading in- 
dicates the region where the gas is molecular, according to Eq. 
1251 Increasingly darker shading shows the regions where carbon 
is found as C"*", C, or CO. The top sequence (a) illustrates the 
effect of decreasing metallicity and dust-to-gas ratio on the dis- 
tribution of C+, CO, and H2. Mostly because of the increase 
in A^H required to attain a given Ay, the CO emitting region is 
pushed further into the clump until, for a fixed cloud size, it disap- 
pears at low enough metallicities. The bottom sequence (b) illus- 
trates the effect of changing the clump size or column density at 
a fixed metallicity (adapted from IBolatto, Jackson Ik lngallslll999l : 
Wolfire, Hollenbach & McKec 201^^ 

At depths where the CO line originates the destruc- 
tion of CO proceeds mainly via photodissociation by 
far-ultraviolet radiation. These dissociations occur via 
line transitions between the ground state and predisso- 
ciated excited electronic states. That is, the majority of 
the transitions to these excited electronic states result 
in dissociation rather than CO in bound, excited lev- 
els. The longest wavelength dissociating transition lies at 
A « 1076 A. Thus the dissociation of CO in the ISM oc- 
curs in the narrow far-ultraviolet band between A ~ 1076 
A and the Lyman limit at A « 912 A. These dissociation 
bands for CO overlap the Lyman bands of H2, and the 
Lyman lines of HI. Thus HI and H2 can shield CO. Be- 
cause CO dissociates in line transitions, similar to H2, it 
can also self-shield at sufficiently large CO column den- 
sities. In addition to gas opacity, dust absorption and 
scattering will attenuate the dissociating radiation field. 

A simple expression illustrates the dependence of the 
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CO dissociation rate on these various process: R — 
X-Ro exp(— 7Ay)0 s^^. The first two factors reflect the 
incident radiation field. Rq is the dissociation rate in 
free space for a radiation field with a particular shape 
and X is a scaling factor that measures the strength of 
the radiation field. The exponential factor reflects shield- 
ing by dust, with Ay the line of sight extinction into the 
cloud. The factor 7 includes both cloud geometry and the 
translation from Ay into extinction at the dissociating 
wavelengths, accounting for absorption and scattering by 
grains in the far- ultraviolet. The parameter accounts 
for shiel ding by HI an d H2 gas and CO self-shielding. 

For a iDraind (|1978l) interstellar radiation field, Rq = 
2.6 X 10-1° s-i (jVisser et all [20091 ): Rq wiU vary for ra- 
diation fields with other shapes. For typical ISM grains 
and penetration into a plane parallel layer, 7 = 3.53 
(Ivan Dishoeck. Jonkheid fc van Hemertl[2006h . The pa- 
rameter O is a monotonically decreasing function of the 
HI, H2, and CO column densities and varies between 1 
(for no shielding) and 0. A fit to the line overlap and 
self-sh ielding function, O, can be found in IVisser et al.l 
(|200a). 

Balancing CO format ion and destruction 

IWolfire. Hollenbach fc McKei (l2010l) find an ex- 
pression for the depth into a cloud at which tqo = 1 for 
the J = 1 — > transition: 



AviRco) - 0.102 In 



3.3 X 10^ 



X 



1 



(24) 



Here n is the gas density and x and Z' are the far- 
ultraviolet field and gas phase abundances in units of 
the local Galactic values. Av{Rco) is the depth in units 
of magnitudes of visual extinction due to dust at which 
Tco = 1- In this expression the explicit dependence on 
dust-to-gas ratio cancels out, but the conversion from Ay 
to column density will depend on the dust-to-gas ratio. 

Variations in Xqq will depend on the relative 
extent of the CO and H2 lay ers. Therefore, 
IWolfire. Hollenbach fc McKed ( 20101) present a similar 
expression for the depth at which the gas is half molec- 
ular 



Av{Rh2) ^ 0.142 In 



5.2 X lO^J^^GR, 



X 



6^ 



DGR 



1.75 



+ 1 



(25) 

where ^dqr is the dust-to-gas ratio relative to the local 
Galactic value. 

The difference between the depth of the H2 and CO 
layers is thus: AAy = 0.53-0.045 In ^-f 0.107 ^JJ^Qp,)- 
0.204 In(Z'). For a dust-to-gas ratio that scales as the gas 
phase metallicity we have 

AAv = 0.53 - 0.045 In - - 0.097 In(Z') . (26) 
n 

The equations for Av{Rco), Av{Rh2), and A Ay all 
show weak dependencies on x, '^dgR' ^^"^ ^' • Thus they 
may be expected to change only weakly with changing 
local conditions and at lower metallicity, the location of 
the Tco = 1 surface, the H2 transition, and the spacing 
between the two all remain approximately fixed in units 
of visual extinction due to du st. The ratio of d ust to 
gas also drop with metallicity (! Draine et al.ll2007l ). As a 
consequence of lower dust abundance and approximately 



constant A Ay , the physical depth of the Tco — 1 surface 
shifts deeper into the cloud, producing a larger surface 
layer of H2 at low Ay- A correspondingly larger layer 
of HI also exists beyond Ay (i?H2 ) but does not bear di- 
rectly on this review. Thus, with lower metallicity and 
correspondingly lower dust abundance, CO retreats fur- 
ther into the cloud than H2 . Figure |5| illustrates the ex- 
pected interplay of HI, H2, CO, and C"*" in a spherical 
cloud and the impact of changing metallicity (top row) 
and cloud size or column density (bottom row). 

If il/(i?H2) is the molecular mass within the radius 
where the molecular fraction is 0.5, and M(CO) is the 
mass within the radius of the tco = 1 surface, then the 
"CO-faint" gas fraction can be defined as / = [M(i?H2 ) ^ 
M{Rco)]/M{Rh2)]- For a cloud with an r'^ density 
profile, corresponding to M(r) c>c r^, the "CO-faint" gas 
mass fraction is given by 



/ = 1 - exp 



-0.76AA 



= — = 1 - exp = 



-A.OAAi 



Ai 



(27) 

where AAy is the optical depth in the "CO-faint" gas 
layer (Eq. |26|) . iV22 is the mean H column density in 



the CO portion of the cloud in units of 10 cm , and 



A 



The retreat of the CO emi tting surface was first noted 
by ea rly studies (|Malonev fc Black 1988; Leq ueux et al.l 
[T99I . This shrinking CO core leads to lower beam fiU- 
ing factors, and thus a lower observed CO iiitensit y (e.g., 
iPak et al.lll998l: iBolatto. Jackson fc Ingallsl [19991) . The 
result is that Xco increases with decreasing metal abun- 
dance. iMalonev fc Blackl (|1988D find nearly equal CO 



and H2 cloud sizes at solar metallicity and so note that in- 
creasing CO abundance above solar will not substantially 
change Xqo- Recent observational estimates find the 
fraction of "CO-faint" molecular gas mass to be ~ 50% 
in the Milky Way (see ii4.2.4p . This allows slightly more 
room for a changing Xco moving to supersolar metallic- 
ities. For example, a change from Xco. 20 ~ 2 to ~ 1.3 
might be plausible due to the effect of increasing metal- 
licity. 

A secondary consideration is that the brightness tem- 
perature of CO may also be affected by systematic exci- 
tation changes with metallicity. The temperature of the 
gas represents a balance between heating and cooling. If 
grain photoelectric heating dominates at the tco = 1 
surface, then to first order both the heating rate (pro- 
portional to the dust-to-gas ratio for photoelectric heat- 
ing) and the cooling rate (proportional to metal abun- 
dance for line cooling) will scale similarly with metal- 
licity, yielding no change in gas temperature with metal 
abundance. It remains unclear, however, whether the 
photoelectric heating efficiency increases, decreases, or 
stays the same at low met allicities (e.z.. [R511ig et al.l 
[200i HsraeT fc Malonevl[20Tll) . 

In summary, models predict an increase in Xco for 
large regions of low metallicity galaxies due to the con- 
traction of the CO-emitting surface relative to the area 
where the gas is H2 for a fixed cloud size. This effect 
may be offset (or compounded) by a mild increase (or 
decrease) in brightness temperature brought about by 
changes in heating sources, chiefly the photoelectric heat- 
ing effect. Note that the usefulness of CO as a tracer of 



24 



total H2 mass will ultimately break down at low enough 
metallicity, where it will be found only at the highest 
column densities and in well-shielded environments. 

6.2. Metallicity-Dependent Calibrations of Xqo 

Calibrating Xco as a function of metallicity has rep- 
resented a key goal of extragalactic CO studies for two 
decades. The issue remains complicated for several rea- 
sons: some practical, some theoretical. We have al- 
ready discussed the biases and caveats of the different 
Xco extragalactic estimators. On practical grounds, the 
faintness of CO in low luminosity, low metallicity sys- 
tems renders observations very difficult. The abscissa 
in any calibration, metallicity, also remains one of the 
hardest quantities in extragalactic astronomy to mea- 
sure with precision. The basic cloud structure arguments 
presented in H6.ll suggest that the sub- and supersolar 
regimes should be treated differently — or at least that 
a single power law represents a poor choice across all 
regimes. We will discuss here several existing calibra- 
tions of Xco with metallicity and contrast them against 
the data, presenting a simple tentative formula that in- 
cludes metallicity effects in SjH 

Some of the most comprehensive and widely 
used work comparing Xco to metallicity relies 
mostly o n vir ial mass-based Xgo estimate s (e.g. 
WilsonI 119951: lA rimoto. Sofue fc Ts uiimotol 11996 ' 



Boselli. Lequeux fc Gavazzi 2002; Bolat to et al.l 12005 



These calibrations should best be viewed as calibrations 
of Xco within the CO- emitting region, and not be used 
to predict the total H2 content from CO on large scales. 
W e suggest that t he lac k of an Xco vs. metallici ty trend 
in iBolatto et al.l ()2008D and the weak trend in (WilsonI 
(|1995l ) reflect the uniformity of opaque, bright CO- 
emitting structures across many types of g alaxies. Sim- 
ilarly the c alibrations bv lArimoto. S ofue fc Tsuiimotol 
(|1996l) and iBoselli. Lequeux fc Gavazzii ()200"^ rely on 
virial mass es to derive Xco as a function of metallicity 
(and so do lObreschkow fc Rawlingil2009D . 

What are the appropriate calibrations to consider 
for Xco as a fun ction o f met allicity? Dust-based ap- 
proaches s uch as llsraell (I1997D or the "alternative ap- 
proach" of IBoselli. Lequeux fc Gavaz zi ( 200 2|) or com- 
bined CO [CII], and dust modeling ([Madden et al.lll997l: 
iPak et aLl 119981) appear to offer the best, or at least 
the most accessible, extragalactic approach. These have 
the potential to capture the whole H2 distribution. In 
the following pagraphs we employ these observations 
to test recent theoretical work on the subject (e.g. , 
Wolfire. HoUenbach fc McKed IMot IGlover et al.l iMa 



Narayanan et al.l 120121 : iFeldmann. Gnedin fc Kravtsovl 
201211 . Note that, as discussed in §6.11 we do not expect 



A'co to change significantly due to metallicity in the su- 
persolar metallicity regime. Instead, variations due to 
optical depth, line width, or excitation temperature are 
likely to dominate. 

Figure [9] presents an attempt at such a compari- 
son. We show an ensem ble of dust and FIR deter- 
mina tions (Israel 1997; Madden et al.l [19971: ILerov et al 



2011 



2007t FGraticr et al. 2010l: lLerov et al.ll2011l:IBolatto et a! 



iSmith et al.. ,20121) target ing local galaxies alon; 



with Xco estimates determined iSandstrom et al.l (|201 
for a sample of 22 nearby disk galaxies. In the 
right panel we plot theoretical predictions for Xco 



as a function of metal abundance, assuming dust-to- 
gas ratio is linearly depend e nt on metallicity. For 
iWolfire. HoUenbach fc McKed (|2010[) we plot the equa- 
tion 



-4.0 AAv -4.0 AAv 
exp ^— exp —= , (28) 
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V.MW 



A 



V,MW 



Xco jZ') 
Xco {Z' = l) 

which is obtained from Eq. [27] assuming that dust-to- 
gas ratio tracks metallicity. Here Z' is the abundance of 
heavy elements and dust relative to solar, Xco {Z' = 1) 
is the CO-to-H2 conversion factor at solar metallicity, 
and Av,MW is the mean extinction through a GMC 
at Milky Way metallicity {Av,mw ~ 5 for Sgmc ~ 
lOOM0pc"2). The prediction for IGlover fc Mac Lowl 
([201 Ih simply adapts their Equation 16, which gives Xco 
as a function of Ay- We assume that the mean ex- 
tinction through a cloud scales as the metal abundance. 

Ay = Z'AvMW 

Xco jZ') ^ r _ 1 for Z'AvMW > 3.5 

Xco {Z' = 1) 1 (Z'Ay.Miy)"' ' for Z'Ay.Mw < 3.5 

(29) 

In the full prediction of iNaravanan et "all ([20121 ). Xco 
depends on a combination of Ico and Z, rather than 
metallicity alone. We cannot readily place their predic- 
tions in this plot. Instead we plot their Equation 6, in 
which Xco depends on Z' and EH2- For each of these 
predictions, we assume SH2 = 100 Mq pc^^, translating 
this t o Ayx) assuming Ry = 3 .1 and the dust-to-gas ra- 
tio of lBohiin. Savage fc Draki ([l9 78). We also plot the 
metallicitv-onlv fit by llsraell ([19971 note that his fit in- 
cludes the radiation field, with a definition that renders 
it heavily degenerate with the dust-to-gas ratio). 

Figure M shows large scatter, even among de- 
terminations using similar techniques for the same 
galaxy, but provides reasonable evidence for an up- 
turn in Xco with decreasing metallicity below Z' ^ 
1/3-1/2. The fit by llsraell (|1997[ ) skirts the up- 
per envelo pe of measured values, whi l e the pre- 
diction by IFeldmann. Gnedin fc Kravtsovl ([2012D pre- 
dicts weaker-than-observed variations at very low 
Z', though data remain scarce (those points rep- 
resent only ~ 4 galaxies). The sharp nonlinear 
increases predicted by shield ing-or iented models like 
Wolfire, HoUenbach fc McKei ([2010f ) or the simulations 
of Gl over fc Mac Low. (,20 111 ) provide the best fits to the 
existing data. 

Less direct approaches to constrain Xco also exist. 
A recently popular technique is to assume an underly- 
ing relationship between star formation and molecular 
gas and to use this relationship and an estimate of the 
recent star formatio n rate to arrive at the H2 present. 
iSchruba et al.l ([2012D applied this technique to estimate 
Xco as function of metalli city in a sample of nearby 
galaxies. iGe nzel et al.l ([20121 ) combined a wide set of low 
and high redshift measurements to estimate the depen- 
dence of Xco on metallicity. Both studies find a signifi- 
cant dependence of Xco on metallicity but with a wide 
range of possible power law exponents, Xco oc Z^ to 
Z^. The strength of this approach is that the observa- 
tions needed to make such estimates are widely acces- 
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Fig. 9. — Conversion factor, estimated from dust-based approaches, as a function of gas-phase abundance. (Left) Color points 
show estimates for v ery n e arby g alaxies Israel (1997), Madden et al. (1997, based on [CII] ), Leroy et al. (2007), IGratier et al.l 0010!), 
IRoman-Duval et al .t (12010!), 'Leroy ct al. (2011), Bolatto ct al. (2011), and Smith et al. (2012). Gray points show high quality solutions 
from analysis of 22 nearby disk galaxies by (Sandstrom et al. 2012), with typical uncertainties illustrated by the error bars near the bot- 
tom left corner. Metallicities are from Israel (1997), Bola tto et al., (2008 ), and M oustakas et al., 12010 ) and quoted relative to solar in 
the relev ant syste m (12 -|- log[0/H] = 8.7 for the first two, 12 -|- log[0/H] = 8.5 for the latter which uses the metallicity calibration by 
IPilvugiii & Thuanl [20051 ). Note that significant systematic uncertainty is associated with the x-axis. The color bands illustrate our recom- 
mended ranges in «co for the Milky Way and ULIRGs. (Right) Colored lines indicate predictions for Xqq as a function of metallicity 
from the references indicated, normalized to Xqq 20 = 2 at solar metallicity where necessary. For these predictions, we assume that GMCs 
have (Sgmc) = 100 M0 pc~^, which we translate to a mean extinction through the cloud using Eq. 1211 Dust-based determinations find 
a sharp increase in Xqq with decreasing metallicity below Z ~ 1/3-1/2 Zq. 



sible. The weakness, of course, is that it requires as- 
suming an underlying relationship between H2 and star 
formation. Any true dependence of the star formation 
efficiency on metallicity, or any other quantity covariant 
with metallicity, will be recast as additional variations in 

If Xqo does increase rapidly moving to low metallic- 
ity, our knowledge of the distribution function of molec- 
ular column density will present a practical limit to the 
usefulness of CO to trace H2. At metallicities perhaps 
as high one half solar, half of the H2 mass will exist 
outside the CO-emitting surface, and that fraction will 
rapidly increase for decreasing metallicity. Thus, appli- 
cation of a CO-to-H2 conversion factor at very low metal- 
licity ultimately involves extrapolating the total mass 
of a cloud from only a small inner part; by Z' ^ 0.1 
this may already be analogous to measuring the total 
H2 mass of a Milky Way cloud from HCN emission or 
some other high density tracer. In the Milky Way, the 
fraction of cloud mass in high column density lines of 
sight can vary dr amati cally from cloud to cloud (e.g . , 
IKainulainen et al.ll2009l: iLada . Lom bardi fc Alv"esll2Q10[ ). 
In the LMC. Ilsrael et al" 111996) found large region-to- 
region variations in the [ClI]-to-CO ratio. As one moves 
to increasingly low metallicities, the use of CO emission 
to quantify the H2 reservoir becomes more and more ex- 
trapolative. While difficult to quantify, this effect should 
add significant scatter that increases with decreasing 
metallicity, and will eventually present a practical floor 
past which CO is not a useful tracer of total H2 mass. 

7. Xco IN STARBURSTS AND OTHER LUMINOUS 
GALAXIES 

Molecular gas in starbursts exists under conditions 
very different from those found in most normal galaxies. 



Observations of starbursts suggest widespread gas vol- 
ume and column densitie s much higher than those typi- 
cal o f normal disks (e.g., IJackson et al.|[T995l : llono et all 
12001 . Molecular gas in starbursts is also warmer, ex- 
citing higher rotationa l transitions than those found in 
less active objects Ce.g..lBradford et al.ll2003l : iWard et"al] 
l2003t iRangwala et al.ll2011f) . In fact, a negative correla- 
tion is observed between molecular gas depletion time 
(a parameter that characterizes how long a galaxy can 
maintain its current star formation rate), and excita- 
tion fe.g.. lPaglione. Jackson fc Ishizukilll997[ ). Similarly, 
a positive correlation exists b etween gas density and star 
formation rate (e.g.. [Gao k, Solomon 2004). These ob- 
servations show that there is a fundamental relation be- 
tween the density and temperature of the molecular gas 
and the existence of starburst activity, such that the gas 
present in starbursting galaxies or regions has higher den- 
sities and temperatures than those prevalent in quiescent 
systems. 

What are the effects of higher temperatures, densities, 
and column densities on Xco? To first order, higher 
gas temperatures yield brighter CO emission, decreasing 
Xqq. Note, however, that while increasing the temper- 
ature decreases A^cOi increasing the density and surface 
density of the self-gravitating clouds of gas has the op- 
posite effect, increasing Xqq (see Eqs. [11] [T^l and [T5|) . 
Therefore we expect a certain level of compensation to 
occur, lessening the impact of environment on the con- 
version factor as long as most of the CO emission arises 
from GMCs. In regions where the average gas density 
is comparable to that of a CMC, however, the entire 
medium will turn molecular and CO emission will origi- 
nate from an extended warm phase. These conditions are 
thought to be prevalent throughout the active regions of 
the brightest starbursts, such as ULIRGs. Adding to 
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this fact, many of these luminous galaxies are mergers, 
or the starbursts occur in regions such as galaxy centers. 
The gas correspondingly experiences motions in excess of 
the velocity dispersion due to its self-gravity. The large 
column densities conspire to make the medium globally 
optically thick, thus setting up the conditions discussed 
in In this situation the CO emission will be dispro- 
portionally luminous, driving Xqo to lower values. 

We can quantitatively explore this scenario 
with PDR model calcul a tions (adapted from 
iWolfire. Hollenbach fc McKed l2010f ). which incor- 
porate self-consistently the chemistry, heating, and 
radiative transfer. Setting up a "typical" Milky Way 
GMC (a virialized structure of size ~ 30 pc with 
Emoi w 170 Mq pc-2 and A/ w 1 x 10^ Mq ) we 
reproduce a Galactic Xqo- In a ^ 160 pc virialized 
cloud with Smoi « lO'' Mq pc-2 and Mw 2 x 10^ Mq 
(representing the molecular structures observed in 
existing high-resolution observations of ULIRGs) under 
the same Galactic conditions, we obtain Xco,20 ^ 80 
in rough agreement with Eq. 1111 Increasing the gas 
velocity dispersion to include 2 x 10^ Mq of stars (see 
and increasing the UV and cosmic ray fluxes 
by 10^ to account for the larger SFR, decreases the 
conversion factor to Xco,20 ~ 0.6. Most of this effect 
is due to the velocity dispersion: increasing the UV 
and cosmic ray fluxes by only 10 with respect to the 
Galactic case yields Xco,20 ~ 2, still much lower than 
t he starting value of 80. 

[Narayanan et al.l (j201lD explore these effects in detail 
using a series of computational models of disks and merg- 
ing galaxies. They find that Xqo drops throughout the 
actively star-forming area in merger-driven starbursts 
due to increased gas temperatures (caused in part by 
collisional thermal coupling between dust and gas, which 
occurs at high densities) , and the very large velocity dis- 
persion in the gas (in excess of self-gravity) which per- 
sists for at least a dynamical time after the burst. The 
magnitude of the drop in Xqo depends on the param- 
eters of the merger, with large Xqo corresponding to 
low peak SFR. Thus, the large drop in Xco occurs in 
massive mergers during the starburst phase, and Xqo 
settles to normal values when the star formation activity 
and the conditions that caused it subside. The simu- 
lated normal disk galaxies experience less extreme con- 
ditions of density and turbulence, and accordingly pos- 
sess higher values of Xqo except in their centers. The 
mean values of the CO-to-Il2 factor for the simulated 
mergers and disks are Xco, 20 ~ 0.6 and 4 respectively, 
with a very br oad distribution f or the mergers. In a 
follow-up study [Narayanan et al.l (|2012f ) introduce a cal- 
ibration of Xco where it becomes a function of W(CO), 
as well as metallicity, Z. This calibration captures the 
fact that the factors that cause a drop in Xqo occur in- 
creasingly at higher surface densities. Note, however, 
that because Xqo is a non-linear function of W(CO) 
obtained from luminosity-weighted simulations, the cal- 
ibration must be applied carefully to observations. The 
observed CO intensity corresponds to the intrinsic in- 
tensity multiplied by a filling factor fbeam < 1, while 
the luminosity-weighted W(CO) employed by the model- 
derived calibration is much closer to the intrinsic CO 
intensity, which is not directly observed because obser- 



vations do not completely resolve the source. 

In the following sections we will discuss the observa- 
tional findings in starburst galaxies in the local universe. 

7.1. Luminous Infrared Galaxies 

One of the earliest comprehensive studies of the state of 
the molecular gas and the value of Xco in a luminous in- 
frared galaxy (LI RG) was pe r forme d on the prototypical 
starburst M 82. IWild et al.l ()1992l) use multi-transition 
observations and detailed excitation calculations to de- 
termine the proportionality between the optically thin 
C^^O emission and the H2 column density. Bootstrap- 
ping from this result, they find that Xco varies across 
the disk of M 82, and is in the range Xco. 20 ~ 0-5 — 1. 
They attribute this low Xco mostly to high tempera- 
tures, as they estimate Tkm > 40 K for the bulk of 
the gas from their modeli n g. Us ing a similar technique, 
iPapadopoulos fc SeaauistI ((1999) analyze the inner re- 
gion of NGC 1068, a well known starburst with a Seyfert 
2 nucleus. They conclude that a diffuse, warm, molecu- 
lar phase dominates the ^^CO emission, while the mass 
is dominated by a denser phase that is better observed 
in the C^^O isotopologue. The ^^CO J = 1 — >■ emission 
from the diffuse phase has low optical depth (ti ~ 1 — 2) 
and is not virialized, and is thus overluminous with re- 
spect to its mass. They argue that the CO-to-H2 conver- 
sion factor in the nuclear region is Xco, 20 ^ 0.2 — 0.4, 
although the precise value depends critically on the as- 
sumed CO abundance; 

Along similar lines, iZhu. Seaauist fc Kunol ()2003D per- 
form a multi-transition excitation study of molecu- 
lar gas in the Antennae pair of interacting galaxies 
(NGC 4038/9). They find that the H2 gas mass in 
their analysis depends critically on the ^^C0/H2 abun- 
dance ratio, which is an input to their model. They con- 
clude that ^^C0/Il2^ 0.5 — 1 X 10^"*, and consequently 
Xco,20 0.2 - 0.4 for the center of NGC 4038 and 
Xco. 20 ^ 0.5 — 1 for the "overlap" region between both 
galaxies. They argue that the adoption of this smaller- 
than-standard Xco is consistent with the gas distribu- 
tion, including HI, in the interacting pair, while the very 
large mass for the "overlap" region that would result 
from adopting a standard Xco would be very difficult 
to reconcile with the gas dynamics. The authors find 
that their excitation analysis points to gas with high ve- 
locity dispersion, large filling fraction, and low optical 
depth as the reas on why CO i s over luminous in the An- 
tennae. Recentlv. iSliwa et al.l (|2012[ ) find also a low Xco 
in their excitation and dynamical analysis of Arp 299 us- 
ing interferometric data, Xco, 20 ~ 0.2 — 0.6. 

The broad conclusions of these in-depth studies are 
in agreement with findings fr om studies carried over 
large samples. lYao et all ()2003D survey 60 local infrared- 
luminous (starburst) galaxies spanning FIR luminosities 
Lfir ~ 10^ - 10^2 Lq, analyzing their ^^CO J = 1 ^ 
and J = 3 — 7- 2 emission. Using the dust temperature as 
a proxy for the gas temperature, and assuming coexten- 
sive emission and a CO/H2 abundance ratio similar to 
what would be expected if all the gas phase carbon was 
locked in CO molecules for Milky Way abundances, they 
conclude that Xco, 20 ~ 0.3 — 0.8 is consistent with the 
available data. 

iPapadopoulos et al.l (|2012D perform a detailed analysis 
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on another large sample of luminous infrared galaxies, 
with ipiR ^ 10^^ Lq. They find that one-phase radia- 
tive transfer models generally match the observations for 
the lower J transitions of CO with a typical Xco,20 ~ 0.3 
{aco ~ 0.6 ± 0.2 Mq (K km s'^ p c^)~M, a value com- 
patib le with previous studies (e.g., iDownes fc SolomonI 
119981 ). The authors conclude that although the gas 
temperature is partially responsible for the lower Xqo, 
the most important factor is the gas velocity disper- 
sion. They also point to a caveat with this result, which 
has to do with the possible existence of a dense, bound 
phase with lower velocity dispersion and a much higher 
Xco- This phase is associated with dense gas tracers 
such as high-J CO or heavy rotor molecules (HCN, CS, 
HCO+). It has the potential to dominate the molec- 
ular mass of the system, and raise the conversion fac- 
tor to Xco, 20 2 — 6. Because the uniqueness of this 
explanation for the observed excitation is unclear, and 
the derived masses conflict in some cases with dynam- 
ical mass estimates, we consider the results of multi- 
compo nent models an inter esting topic for further re- 
search. iPapadopoulos et al.l also point out that although 
the global CO J = 1 — >■ luminosity is dominated by 
the warm, low Xqo component, it is possible to hide 
a cold, normal Xqo component that could add a signif- 
icant contribution to the molecular mass of the system. 
This component would most likely be spatially distinct, 
for example an extended molecular disk, and thus could 
be separated in high spatial resolution studies with in- 
terferometers. 

7.2. Ultraluminous Infrared Galaxies 

Ultra-luminous infrared galaxies (ULIRGs) are ex- 
treme cases of dust-enshrouded stabursts and active 
galactic nuclei, with far-infrared luminosities in excess 
of LpiR ~ 10^^ Lq. These objects are the products of 
gas-rich m ajor mergers, and posses s very large CO lu- 
mino sities (jMirabel fc Sanderj [19881 iSanders fc Mirabell 
119961) . Despite that, their Lfir/Lco ratios frequently ex- 
ceed those found in spiral galaxies, including interacting 
pairs, and Milky Way GMCs. The fundamental question 
is how those CO luminosities relate to their molecular gas 
masses. Understanding Xco in ULIRGs is particularly 
interesting because they likely provide the best local tem- 
plates for the most luminous high-redshift submillimeter 
galaxies, which have very high star formation rates com- 
pared to other galaxies of similar stellar mass at their 
redshift (i.e., are off the m ain sequence, iTacconi et al.l 
[200l [Narayanan et al.|[2010f ). 

Studies presented in a ser ies of papers 
jDownes, Solomon & Radford' '1993; 'Solomon et al.l 
1997; Downcs & Solomon 1998; Brvant & Scovillc 1996, 
1999) show that the molecular gas masses obtained for 
ULIRGs using the Galactic disk CO-to-H2 conversion 
factor are uncomfortably close to (or exceed) their 
dynamical masses, suggesting that the Galactic Xco 
overestimates their total molecular mass. These authors 
develop a consistent one-component model that explains 
the high-resolution observations as rotating, highly 
turbulent remnants of the merging process. The CO 
emission is dominated by low density gas, and although 
it is opt ically thick it is o nly moderately so (for example, 
see also llono et all I2007D . This lower molecular mass is 
consistent with the observed (optically thin) millime- 



ter dust continuum emission for Galactic dust-to-gas 
ratio. The "typical" Xco they derive is Xco, 20 ^ 0.4, 
approximately a factor of 5 lower than in the Milky 
Way disk (aco ^ 0.8 M© (K km s'^ pc^)^!, with 
individual results ranging between 0.3 and 1.3 in Table 
9; .Downes fc Solomon 1998). 

Detailed high res olution studies of individual U LIRGs 
find similar results. iBrvant fc Scovilld (|1996lll999[ ) study 
the kinematics of seven LIRGs and ULIRGs, determining 
interferometric sizes and dynamical masses. They con- 
clude that the use of a Milky Way conversion factor re- 
sults in molecular masses larger than the dynamical mass 
of the system in the cores of all seven objects, although in 
some cases that can be explained away assuming a face on 
orientation. They find Xco, 20 < 0.7 and Xco, 20 < 1-5 
for Mrk 231 and NGC 6240 respectively Modeling 
the CO kinematics of Arp 220, including the effects of 
the s tellar components in th e dyna mics of the nuclear 
disk, iScoville. Yun fc BrvantI (|1997[ ) find Xco,20 ~ 1- 
These estimates are comparable with those from one- 
component radiative models for the same galaxies, and 
typically much lower than the two-component results 
(Downes fc S olomon 1998; Papadopoulos et al., ,201^) . 

7.3. Synthesis: Xco «^ Starbursts 

There has been much recent progress, through both ob- 
servations and modeling, on the determination of Xco 
in starbursts. Studies of CO excitation in large samples, 
as well as detailed studies of individual cases, strongly 
point to lower than Galactic values of Xcoi particu- 
larly in extreme starbursts such as ULIRGs. These 
low values of Xco are driven by a globally molecu- 
lar medium coupled with high gas temperatures and, 
more importantly, very high velocity dispersions in the 
CO emitting gas due to a combin ation of merger ac- 
tivity and the stellar pote ntial (N aravanan et al] 120111 ; 
IPapadopoulos et al.l I2012D . The standard practice 
has been to adopt aco y 0.8 Mq (K km s^^ pc'^)"^ 
(jDownes fc SolomonI Il998f) . a value similar to the 
aco « 0.6 ± 0.2 Mq (K km s'^ pc^)-^ resulting from 

recent one-compone nt modeling in a large sample 

([Papadopoulos et al.l I2012D . Note, however, that this 
is an average val ue and there are likely large galaxy-to- 
galaxy variations (jNaravanan et aT1l201lD . To first order 
these variations should be mostly related to the total sur- 
face density of the galaxy, if the gas is bound and expe- 
riencing the global gravitational potential (c.f., Eq. [T6\\ . 
We will use this to suggest a tentative Xco correction in 

m 

Large uncertainties exist in the estimates from observa- 
tions, stemming from assumptions about ^^C0/H2 and 
^^CO/^'^CO ratios, coextensive emission in the different 
transitions, and the general need for simplifications in 
excitation models. Note also, as a persistent caveat, 
the possibility of hiding significant molecular mass in 
ULIRGs in a low velocity dispersion, dense component. 
This posibility, however, seems disfavored in some cases 
where high quality dynamical estimates are available. 

8. Xco AT HIGH REDSHIFTS 

The improved sensitivity of millimeter and centimeter- 
wave instruments has allowed detection of CO emission 
from an increasingly diverse range of systems out to 
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higher and higher redshifts. Most of the objects ob- 
served in CO to date represent the bright, rare end 
of the luminosity distribution: so-called submillimeter 
galaxies (SMGs) and QSO hosts. Some recent studies, 
however, target rotating disks that lie at the high-mass 
end of the star-forming blue sequence at their redshift 
(i.e., have typical star formation rates given their stel- 
lar mass — we will refer to them as "main sequence" 
galaxies) . CO observations thus now begin to sample the 
regime of "inain s equence" galaxies (|Tacconi et al.ll20Tol : 
iDaddi et al.l 120101 ) . and will expand to lower luminosity 
systems over the next decade. 

Estimating molecular gas masses is frequently the main 
goal of CO observations at high redshifts. Unfortunately, 
direct determination of Xqq in high redshift objects re- 
mains tremendously challenging. An additional compli- 
cating factor is that most high-redshift observations do 
not measure the J = 1 — ^ transition of CO, but higher 
rotational transitions. Thus translating these measure- 
ments into molecular masses requires, at least in princi- 
ple, understanding CO excitation in addition to the other 
physics governing Xco- 

At this st age, the modeling of opti cally thin isotopo- 
logues (e.g.. lPapadopoulos et al.|[2012D may offer the best 
opportunity for direct Xco measurements at high red- 
shift. Such observations have been too costly to under- 
take, and even with ALMA may only be practical for 
the brightest objects. Dust continuum observations of- 
fer another route to estimate molecular masses, although 
their application at high redshift requires understanding 
dust-to-gas ratios in relation to metallicities — a complex 
problem that ultimately requires knowledge of the bal- 
ance of ISM enrichment, dust production, and destruc- 
tion processes in galaxies. In lieu of direct measurements, 
the best route is to understand the physical drivers of 
Xco and to apply knowledge acquired from local galax- 
ies to systems at high redshift. 

8.1. Observed CO Line Ratios 

The ratio of the J = 3— >2toJ=l— >0 tran- 
sition, rai = T3/T1, exhibits a wide range of values 
across galaxie s and has been obse r ved a t bot h low and 
high r edshift. iMauersberger et all (|1999[ ) and lYao et ahl 
(|2003[ ) consider r^i for large samples of nearby galaxies 
with infrared luminosities 9 < log(LFiR) ^ 12. They 
find a mean r^i w 0.63-0.66 (median rai 0.5) and 
a broad distribution of values. No strong correlations 
between rai and dust temperature or luminosity are ev- 
ident, but they note that rai increases with increasing 
concentration and star formation efficiency (measured 
as -^ftr /M(H2 ), though n ote the CO intensity in the 
denominator). Ilono et all ((2009) examine 15 luminous 
LIRGs and ULIRGs with \og{Lpiji) > 11.5 and find mean 
r^i « 0.48 (median 0.4) integrated over whole galax- 
ies, with a higher mean r^i ft- 0.96 a t the location of 
peak CO emission. iMao et all ()2010f ) discuss 7-31 in a 
sample of over 60 galaxies, finding that barred galax- 
ies and starbursts have the higher averaged r^i values 
(rai ~ 0.89 ± 0.11 for starbursts) followed by AGNs 
see t heir Table 4). More recently, IPapadopoulos et all 
20 111) report the results of a large sample of LIRGs 
and ULIRGs composed of new observations and liter- 
ature compilation, where they find mean r^i ss 0.67 and 
r2i « 0.91. By contrast with these samples of active 



galaxies, the inner po rtion of the Milky Way has a ra - 
tio rai w 0.28 ± 0.17 (jFixsen. Bennett fc MathCTlll999D . 
denoting a correspondingly lower excitation in a z = 
"normal" star forming galaxy. 

To first order r^i offers a tool to distinguish excited, 
star-forming gas similar to that seen in LIRGS and 
ULIRGS from the gas found in more normal systems 
like the Milky Way. As emphasized by a number of 
authors, however, these line ratios are hardly unam- 
biguous indicators o f loca l conditions. For example, 
IPapadopoulos et ahl (j201lD show that for their sample 
of LIRGs and ULIRGs possible densities and temper- 
atures are often poorly constrained, with data allow- 
ing high density {n ^ 10^ cm"'^), normal temperature 
{Tkin ~ 20 K) or low density {n ^ 300 cm~^) high tem- 
perature (Tkin ~ 150 K) solutions. Note as a limitation 
the assumption of similar filling fractions for galaxy-wide 
averages in all transitions, while high-resolution maps 
frequently find excitation gradients. It is thus key to 
appreciate that r^i, although useful as a first order exci- 
tation indicator, combines the effects of both excitation 
and beam filling. 

With the improving sensitivity of the available in- 
strumentation, measurements of line ratios including 
CO J = 1 — >■ at high redshift are rapidly increasing. 
For SMGs, the presumed high redshift simile of lo- 
cal ULIRGs, the measured line ratios resemble those 
found for local LIRGs and U LIRGs. Harr is ct al. (2010), 
iSwinbank et al.l (|2010D . and llvison et al.l (|2011[ ) report 
mean r^i « 0.6 8, 0.66, and . 55 res pecti vely for a hand- 
ful of SMGs. ICarilh et al.l (IMol) and IBothwell et all 
(|2012l ) note that fitting observed SMG spectral line en- 
ergy distributions requires a cool and a wa rm component, 
with most of the mass in the former. iRiechers et al.l 
(2011b) find low integrated r4i « 0.25 — 0.6 in two 
SMGs, noting the existence of excitation gradients, with 
CO J = 1 — >■ in a more extended distribution than the 
higher tr ansition. At the brigh test end of the galaxy dis- 
tribution EiechersiraD (|2011l l report rgi - 0.96 - 1.06 
for four high-redshift SMGs hosting quasars, suggesting 
that the excitation in these sources is such that the emis- 
sion is thermalized to higher rotational transitions (see 
also lWeiss et al.l [20071 for a CO excitation discussion of 
a number of bright SMG and QSO-host high-redshift 
galaxies) . 

The lowest galaxy luminosities currently probed by 
high redshift CO observations correspond to the high 
luminosity end of the star-forming galaxy main se- 
quence. Unlike SMGs or local ULIRGS, these systems 
are not major mergers. Excitation da ta remain scarce for 
these "main sequence" disk galaxies. iDannerbauer et aD 
(2009) observed a BzK galaxy at z ~ 1. 5 report r32 ~ 
0.5 and r2i ^ 0.6. lAravena et al.l (|2010D report r^i « 
0.61 and r2i ~ 0.4—1.2 in a few BzK galaxies. These ob- 
servations have large uncertainties, but suggest that the 
conditions in high-z disk galaxies are not as extre me at 
those in local ULIRGs. As a counterpoint, Ricchcrs et al.l 
([2010^ discuss the CO J = 1 detection of two highly 
magnified (/i ~ 30) low mass galaxies where they find 
7-31 « 0.72 and 0.78 with errors of ~ 0.2, suggesting 
higher excitation than the previous examples if there is 
no differential lensing. For further discussion see the re- 
view by Carilli & Walter in this same issue. 
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8.2. Estimates of Xqo *^ High Redshift Systems 

A number of observational studies have considered the 
problem of Xqo in high-rcdshift sources, combining line 
ratio measurements, consistency arguments, and scaling 
r elations to argue for p lausible values. 

iTacconi et all (|2008f) discuss the values of Xco appli- 
cable to a sample of high-z star-forming galaxies domi- 
nated by SMGs but also including lower mass galaxies. 
Allowing for a dark matter contribution of 10% — 20% 
within a 10 kpc radius, they minimize the of the dif- 
ference between the dynamical mass and the sum of the 
stellar, gas, and dark matter. They find a Galactic Xco 
to be strongly disfavored, while a smaller, ULIRG-like 
-^CO 20 ~ 0.5 produces much more satisfactory results 
for SMGs (Fig. M). 
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Fig. 10. — Minimization of the difference between dynamical 
mass and the sum of the gaseous, stellar, and dark matter com- 
ponents, for a sample of 9 high-z galaxie s with high quality data, 
dominated by SMGs (Tacc oni et ai.l200a) . ^ind reduced (left 
and right side axes, respectively) are shown here as a function of 
chosen Initial Mass Function, for different gas (a) and dark matter 
ifdark) parameters. The analysis distinguishes between a in the 
SMG subsample (asMc) ^^'^ iri the UV selected subsample {«uv) 
of lower mass galaxies. 



In their SMG sample, llvison et all 1)20111 ) find 
-^^00,20 ^ 0.9 — 2.3 to be compatible with the ob- 
served star formation histories and dynamical and stel- 
lar masses. They argue for lower Xco, 20 ^ 0.5 in most 
cases, however, based on star formation efficiency con- 
siderations. They decompose the molecular gas into a 
warm, star-forming phase (with intrinsic rai ^ 0.9) and a 
cold, quiescent phase {r^i ~ 0.3), and invoke an approx- 
imately maximal starburst for the star-forming phase. 
The maximal starburst, Lir/M(H2) « 500 Lq Mq\ 
represents the largest allowed ratio of star formation to 
H2 before radiation pres sure disperses a starburst (e.g., 
iThompsqn. Quataert &: Murray, 2005.'). 

While ITacconi et al.l (|2008D imd" lIvison et al.l (poTTI ) 
poin t to similarities between ULIRGs and SMGs, note 
that iBothwell etall (|2010D find that several SMGs have 
gas distributions (imaged in J > 3 transitions) more 
extended than those of local ULIRGs. This may sug- 



gest that local ULIRGs may not be the best analogues 
of SMGs, and that some SMGs may not aris e from major 
mergers. The excitation g radients seen by llvison et al.l 
( 20li |) and lRiechers et alT (|201l) also point to extended 
quiescent gas reservoirs in SMGs that are not frequently 
o bserved in local ULIR Gs. 

iDaddi et al.l H2010b[ ) and iGenzel et all (|20To[ ) study 
high-redshift galaxies selected from optically identified 
objects at z ^ 1.2 — 2.2. These are still massive sys- 
tems, but drawn from the main sequence at their red- 
shifts rather than major merger-driven starbursts. Kine- 
matically, most of these systems are extended rotat- 
ing disks, alt hough with la r ger v el ocity dispersions tha. n 
loca l disks. IDaddi et al.l (|2010D . IGenzel et ai] ()2010D . 
and ITacconi et al.'r ( 2010D argue that, because their CO 
emission likely arises from collections of self-gravitating 
GMCs (as suggest ed by local estimates of Toomre Q, 
IGenzel et al.ll20l"TI ) with modest temperatures T ^ 20 — 
35 K, these systems have Galactic Xqo- Additional con- 
straints come from dynamical measurements that, to- 
gether with simulations and modeling of the stellar popu- 
lations, yield a gas mass resulting in Xco, 20 ~ 1.7 ± 0.4, 
a value considerably hi gher than observed in ULIRGs 
or estimated for SMGs (jPaddi et al.ll2010b( ). Note that 
these authors also provide a calibration for the dynamical 
mass estimation in this type of galaxy. 

A different approach to quantify the amount of gas in 
a galaxy is to try t o use the continuum dust emission. 
IMagdis et al.l ()2011[ ) use this approach in two very well 
studied high-z galaxies, a massive SMG at z ^ 4 and 
a "main sequence" disk star-forming galaxy at z ^ 1.5. 
They take advantage of the well characterized SEDs of 
these objects to calculate their dust masses employing 
different approaches to dust modeling (traditional mod- 
ified blackbody as well as the type of models used by 
IDraine et al.l I2007f ) . This dust mass is converted into 
a gas mass (assumed to be dominated by molecular 
gas) using a metallicity derived through the correspond- 
ing S FR-mass-metallicity relation fe.g.. iMannucci et al.l 
[2OIOI) and assuming a Galactic dust-to-gas ratio cor- 
rected (approximately linearly) by metallicity. The result 
is A:co,20 ~ 0.4±0.2 for the SMG, and Xco .2n 1.9±1.4 
for the disk galaxy. iMagnelli et al.l ()2012|) follow a sim- 
ilar procedure in a mixed sample of 17 high-z galaxies, 
including "merger-like" and "main sequence" disks, find- 
ing values of Xqo consistent with Galactic for the "main 
sequence" disks and and a factor of several lower for the 
"merger-like" galaxies. They also find Xco « '^dust' 
similar to the behavior expected from Eq. [T^l 

Although this new use of dust observations at high- 
redshift is promising, we caution that going from dust 
SEDs to gas masses involves a large chain of assump- 
tions that is fraught with potential problems. At a fun- 
damental level, we are far from understanding the dust 
creation-destruction b alance which sets the dust-to-gas 
ratio in galaxies (e.g., iDraind [20091 iDwek fc Cherchnefil 
I2OIII and references therein). Moreover, the transla- 
tion of a SED into a dust mass relies almost entirely 
on dust grain properties (composition, mass emissivity, 
size d istribution) that have bee n derived for the Milky 
Way (|Draine fc Lil 120011 l2007bl ). a nd are extremely un - 
certain even for local galaxies fe.g.. lGalliano et"aLll201lD . 
Note that the self-consistent dust-based Xco estimates 
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discussed ! j5.2l are emphasized precisely because they are 
free to a large degree from these problems. Nonetheless, 
this is an interesting approach that reinforces what seems 
to be the outstanding trend discussed above. ULIRG-like 
conversion factors for SMG (merger) galaxies, and Milky 
Way- like conversion factors for disks at high redshift. 
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Fig. 11. — Mctallicity effects in high redshift "main sequence" 
galaxies (Gcnzel et al. 2012). The top panels show rest-frame H« 
images of zC406690 and Q2343-BX610, two 2 ~ 2 galaxies with 
similar rotational velocities (vrot ~ 224 and 216 km respec- 
tively), star formation rates (SFRr^ 480 and 212 Mq yr~^, in- 
cluding extinction correction), and stellar masses (M« ^ 4.3 and 
17 X 10^" Mq). The bottom panels show the corresponding rest 
frame optical spectrum containing the Ha and [Nil] transitions at 
6563 and 6585 A (red), and the CO J = 3 — )■ 2 emission (blue). 
The galaxy zC406690 has a low mctallicity, indicated by a low 
[Nll]/Ha ratio, and a correspondingly low CO emission, despite 
its large Ha flux and star formation activity. 

IGenzel et al.l (|2012[ ) present 44 low-mass high-redshift 
galaxies, showing that observations are starting to probe 
the low mctallicity regime at high redshift. The authors 
find the effect already discussed for local low mctallicity 
galaxies in ^ where CO is disproportionally faint for the 
star formation activity (Fig. [TT]) . By assuming that these 
galaxies obey the same gas surface-density to star for- 
mation relation obse rved in local disks (e.g.. iLerov et al.l 
|2008[) . IGenzel et al.l estimate H2 masses and derive a 
metallicity-dependent calibration for Xco at high red- 
shift. They find Xco oc Z~^-^ — Z~^ -^, in approximate 
agreement with loc al measurements (jLerov et al.ll201ll : 
ISchruba et alllMl) . 

8.3. Synthesis: Xqo o,t High Redshift 

Estimates of the appropriate Xqo for high redshift 
galaxies rely largely on scaling and consistency argu- 
ments. Because they probe the physical state of the 
molecular gas, line ratio measurements represent a par- 
ticularly powerful tool for such arguments. Given the 
difficulties in observing the J = 1 — !• line, they also 
represent a critical measurement to understand how 
to translate observations of high- J transitions into CO 



J = 1 ^ 0. Studies of local LIRGs and ULIRGs show 
^31 ^ 0-5 globally in these sources, and r^i ~ 1 when 
focusing on the active regions. Similar ratios are ob- 
served in very massive galaxies at high redshifts (SMGs), 
where r^i ~ 0.5 — 0.7 for integrated fluxes, and r^i ^ 1 
is observed in the very compact starbursts surround- 
ing QSOs. CO J = 1 — ?> data for main sequence star- 
forming galaxies at high-z remain scarce, making it is dif- 
ficult to assert whether they possess large low excitation 
reservoirs of molecular gas. As the data quality, resolu- 
tion, and availability improves, we expect that the sim- 
ple single-component models used today will evolve into 
more realistic multi-component models for the molecular 
ISM of these sources. 

A number of authors have explored the value of Xco 
applicable to samples of high-redshift galaxies. The con- 
sensu s, validated to first ord er by numerical modeling 
(e.g., INaravanan et al.l 120111) . is that massive merger- 
driven starbursts such as SMGs are most consistent with 
a low Xco similar to local ULIRGs, while blue-sequence 
galaxy disks most likely have higher Xco, similar to lo- 
cal disks. This is reasonable in terms of the physics that 
drive the value of Xco- High density environments with 
an extended warm molecular phase not contained in self- 
gravitating clouds will result in low Xco: while molecular 
gas contained in collections of self-gravitating GMCs will 
have Xco close to the Galactic disk value. We expect the 
most heavily star-forming of these disks to have higher 
II2 temperatures, and consequently somewhat lower Xco 
values. Because of the opposing effects of density and 
temperature in self-gravitating GMCs (for example, Eq. 
[T^ . however, it appears that values as low as those ob- 
served in ULIRGs will only occur when the CO emission 
is dominated by an extended warm component that is 
not self gravitating (e.g., Papadopoulos ct al. 2012). 

Multi-line studies including both low- and high- J tran- 
sitions allow more rigorous constr aints on the density 
and k inetic temperature of the gas (jPapadopoulos et all 
I2OI2D . In the near future, multi-line studies of CO iso- 
topologues or paired observations of CO and optically 
thin sub-mm dust emission will offer additional, more di- 
rect constraints on Xcoi though with their caveats and 
shortcomings. 

An exciting development at high redshift is the emer- 
gence of main sequence "normal" galaxy surveys that 
begin to sample the lower mctallicity regime prevalent 
at early times. In these samples we see Xco effects 
that are consistent with those already discussed for local 
galaxies (®: metal-poor galaxies a re underluminous i n 
CO for their star formation activity (jCenzel et al.ll201^ . 
This situation will be increasingly common, and more 
extreme, as CO surveys sample lower galaxy masses and 
earlier times: by z ~ 3 — 4 galaxy metallicities drop by 
0.5 — 0.7 dex for a fixed stellar mass, and the charac- 
terist ic galaxy mass also becomes lower ()Mannucci et al.l 
|2009() . Ultimately the best probe of the molecular gas 
content of these low mctallicity young galaxies may be 
the fine structure line of [CII], or perhaps the dust con- 
tinuum. 

9. CONCLUSIONS AND OPEN PROBLEMS 

We have summarized the efforts to measure Xco in the 
Milky Way and other galaxies, as well as the theoretical 
arguments and studies that show that CO can be used 
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as a tracer of molecular mass, under certain conditions, 
through the adoption of a CO-to-H2 conversion factor. 

In the following paragraphs, we would like to offer con- 
cise answers to a few key questions that have been de- 
veloped elsewhere in this review: 

• What is the "best" value for the Milky Way Xco? 

• Under what circumstances is the assumption of 
"Galactic Xco" likely to be good or bad? 

• How well can we calibrate Xco as a function of 
metallicity? 

• What is the most appropriate value of Xco in 
ULIRGs? 

• What can we say about Xco in the distant uni- 
verse? 

• Can we offer a practical prescription for Xco based 
on observables? 

• What observations and calculations are necessary 
to move forward? 

Measurements in the Milky Way have achieved a very 
good level of sophistication and consistency, beyond 
what is possible in external galaxies. We have shown 
that there is an assuring degree of uniformity among 
the large scale values of Xco obtained through differ- 
ent techniques, particularly in the inner disk (1 kpc< 
^5 9 kpc). We recommend adoption of a con- 
stant Xco = 2 X 10^° cm-2(K km s^^)"! (aco = 
4.3 Mq (K km s~^ pc^)~^) with an uncertainty of ±0.1 
dex ( a factor of 1.3) for the inner disk of the Milky Way. 
The evidence for a large scale Galactic Xco gradient 
and its magnitude are, presently, at best unclear, and 
the simplicity of a constant conversion factor is prefer- 
able. Nonetheless, there is convincing evidence that Xco 
in the Milky Way center region is smaller than in the 
disk by factors of 3 — 10, and this is reaffirmed by ex- 
tragalactic observations that find that a low Xco is not 
uncommon in other galaxy centers. Following the results 
obtained by a number of studies, we recommend using 
Xco ~ 0.5 X 10^0 cm-2(K km s'^)"^ for R < 500 pc 
in the Milky Way. The uncertainty is difficult to quan- 
tify, but taking the different measurements obtained as 
well as the results discussed for galaxy centers and star- 
bursts, we recommend ±0.3 dex (a factor of 2). We ex- 
pect that further CO multitransition modeling, 7-ray, 
and dust continuum studies will help constrain better 
this value as well as that of a possible large scale Galac- 
tic gradient. For the outer Milky Way we expect Xco to 
increase, in principle, following the same physics under- 
lying its increase in low metallicity environments, as we 
discuss below. 

It is important to recognize that these are average num- 
bers, strictly valid for GMCs on scales of tens of pc. 
The validity of invoking anything like a constant Xco 
on a line-of-sight by line-of-sight basis is considerably 
less defined. Both theoretically and observationally, we 
have shown that considerable variation can exist on small 
scales, reflecting local chemistry and physical conditions. 

The most mature techniques in "normal" galaxies re- 
main virial mass measurements and the use of dust as an 
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Fig. 12. — Conversion factor as a function of total surface density 
for nearby disk galaxies and ULIRGs. The gray points illustrate 
the high S/N solutions for oco based on dust emission on kpc 
scales in a sample of nearby disks, with typical errors illustrated in 
the lower left corner (Sandstrom et al. 2012). The corresponding 
surface densities are dominated by the stellar component. For the 
ULIRGs we plot the aco deter minations by Downes & Solomo^ 
(1998. magenta) and Papadopoulos et ai.l 1 12012 . multitransition 
one-component fits in black and two-component fits in blue, error 
bars represent possible range). The abcisa i s from the dynamical 
mass measurements by 'Downes fc So lomo n| l|1998l ) , thus we plot 
only the overlap of both samples. The color bands illustrate the 
recommended ranges for Milky Way and ULIRG conversion factors. 
The S""'^ line for S > 100 Mq pc~^ is a reasonable representa- 
tion of the overall trends when considering the one-component fits. 



optically thin tracer. Both techniques have their draw- 
backs, and we particularly emphasize the ambiguous (at 
best) sensitivity of virial mass measurements to any ex- 
tended envelope of H2 mixed primarily with C"*" rather 
than CO. At solar metallicities, a wide range of measure- 
ments yield Xco, 20 ~ 1-4, but with large (stiU > factor 
of 2) scatter and uncertainties related to the dynami- 
cal state of clouds and the environmental dependence of 
dust properties. In the absence of further characteriza- 
tion or studies, we recommend the conservative approach 
of adopting Xco = 2 x 10^° cm^^(K km s^^)^^ with an 
uncertainty o/±0.3 dex (a factor of 2) in the disks of nor- 
mal, solar metallicity galaxies. This applies to galaxies 
where the CO emission is dominated by self-gravitating 
II2 clouds or cloud complexes. This value can approxi- 
mately be applied down to metallicities '-^ 0.5Zq, and in 
regions where the total gas plus stars surface density is 
< 300 M0 pc"2. ALMA wih greatly expand the appli- 
cation of both dust and virial mass techniques. Other 
future prospects will include modeling of resolved [CII] 
emission, now widely available thanks to Herschel, the 
extension of spectral line modeling beyond bright galaxy 
centers, and further exploitation of galaxy scaling rela- 
tions. 

Several regimes exhibit clear departures from a Galac- 
tic Xco- Some, but not all galaxy centers share a value 
of Xco lower than disks, qualitatively similar to that 
observed in overwhelmingly molecular massive starburst 
galaxies such as ULIRGs. Dust and spectral line mod- 
eling of these central regions show depressed Xco com- 
pared to galaxy disks, with the depression spanning a 
wide range of Xco up to an order of magnitude below 
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Galactic. These central Xqo depressions are not univer- 
sal, although th ey seem correlated with the stellar sur- 
face density ()Sandstroni et alll2012D . and additional 
information appears necessary to predict the appropriate 
Xqq for use in any specific galaxy center. 

A wide variety of evidence points to high Xqo in low 
metallicity regions: the absolute and normalized faint- 
ness of CO, high [ClI]-to-CO ratios, high SFR-to-CO 
ratios, large dust-based Xco determinations, and the- 
oretical calculations of cloud structure. Virial masses 
represent a significant exception to this body of evi- 
dence. When derived at high spatial resolution, these 
tend to show little or no enhancement in Xqo above the 
Galactic value, even in low metallicity systems. We re- 
view the theoretical expectations for the shrinking of the 
opaque CO emitting surface relative to H2 as metallic- 
ity decreases. We favor a self-consistent use of dust, an 
optically thin tracer of g the currently most ma- 

ture methodology to robustly estimate molecular mass 
at low metallicity. We also highlight the problems with 
the blind use of a dust-to-gas to metallicity calibration. 
The present self-consistent dust-based Xqo estimates of- 
fer support for the picture of rapidly increasing Xco 
at low metallicity, but still yield a wide range of Xqo 
even for similar techniques applied to the same galaxy. 
We recommend adopting a shieldinq-based prescription 
such a s t hat presented bv\Wo lfire. Hollenbach & McKed 
1(20 !(}) or \Glover & Mac Lou (201 li) to account for the 
effects of metallicity, with the normalization chosen to 
match a "Galactic" Xco o,t solar metallicity. 

The uncertainties involved in any metallicity- 
dependent correction remain substantial. As a first 
order picture, we expect Xco will change slowly for 
metallicities larger than 12 -I- log[0/H] ^ 8.4 (ap- 
proximately Zq/2), and considerably faster at lower 
metallicities. By the time Xqo reaches 10 times the 
Galactic value the CO emitting surface encompasses 
only ~ 5-10% of the H2 cloud, suggesting that the 
utility of CO as a global tracer of H2 will become 
more and more marginal as one moves to progressively 
metal-poorer environments. Rather, CO will in fact be 
a tracer of high column density peaks and well-shielded 
regions. 

There is general agreement that the processes oper- 
ating in overwhelmingly molecular and turbulent star- 
bursts, where high gas temperatures are also present, 
drive Xqq to values that can be substantially lower than 
in solar metallicity galactic disks. Because of the role 
of velocity dispersion in setting A"cO: in the absence 
of a self-regulating mechanism it is almost certain that 
there exists a continuum of values rather than a unique 
Xco value that is applicable in all cases. The typ- 
ical result of the one-component modeling is Xqq ^ 
0.4 X 10^° cm~^(K km s~^)~^ for massive, turbulent, ul- 
traluminous starbursts. The range around this value is 
large, at least ±0.5 dex, most of it likely representing real 
physical variation among sources. 

The picture of Xqo at high-redshift is still emerging, 
and instruments like ALMA will make a crucial con- 
tribution to better understand it. Lacking direct mea- 
surements, the best approach is to use knowledge of the 
physical drivers of Xco developed in local galaxies, as 
well as scaling and consistency arguments. The simplest 
approach is to identify the brightest, off-main sequence 



massive SMGs likely due to starbursting mergers with lo- 
cal ULIRGs, while disky, rotation-dominated "main se- 
quence" galaxies are to first order more likely similar 
to local disks dominated by self-gravitating or virialized 
molecular clouds. This is an area of active research. The 
picture will become more nuanced as new observations 
revealing the resolved kinematics of the molecular gas 
and its excitation are obtained. In particular, observa- 
tions of "main sequence" galaxies at 2; ^ 1 — 2 suggests 
that metallicity effects will become an increasingly im- 
portant consideration at high-z, as observations push to 
lower galaxy masses at higher re dshifts and consequ ently 
more metal-poor environments (jCenzel et al.|[20T2l ). 

9.1. Toward a Single Prescription 

Ultimately, we desire a prediction for Xco based 
on observable properties, for objects ranging from 
low-metallicity dwarf galaxies to high surface-density 
ULIRGs. In the following paragraphs we present some 
steps in that direction, referring to aco since that is the 
quantity most often useful for distant galaxies. 

Based on the discussions in and SjTl aco can be 
thought of as having two primary dependencies; one re- 
lated to the temperature and velocity dispersion effects 
driving a low value in ULIRGs, the other related to the 
dominance of CO-faint molecular gas driving a high value 
at low metallicities. Treating the two effects as separa- 
ble, aco = aco,Mw/cOF/sB, where aco.MW represents 
an overall normalization under Milky Way disk condi- 
tions. The factor /cof corresponds to a correction that 
accounts for the fraction of II2 mass associated with the 
outer layers of clouds where most CO is photodissoci- 
ated. The factor /sb accounts for changes in aco due to 
temperature and velocity dispersion. 

Drawing from ^ 36. 11 /cof may be approximated by con- 
sidering Eq. [27] applied to a population of identical, fixed 
surface density clouds, 

/cof 0.67 cxp ( "L^i'oo ) • (30) 

^GMC/ 

Here we assume that dust-to-gas ratio tracks metallicity, 
Z' is the metallicity normalized to the solar value, and 
^GMC i^ the average surface density of molecular clouds 
in units of 100 Mq pc"^. 

The factor fsB is considerably more tentative. The 
simple theoretical arguments we out line in Sj^l as well as 
simulations fe.g.. iShettv et al.]|2011| ). suggest that both 
the gas velocity dispersion and temperature are key pa- 
rameters. Nonetheless, keep in mind that the fundamen- 
tal driver of A"co is what fraction of the CO luminosity 
arises from gas in self-gravitating clouds, versus an ex- 
tended not self-gravitating component bound by the to- 
tal mass of the system. Given current observational con- 
straints and our desire to parametrize in terms of measur- 
able quantities, we suggest that the variations between 
normal disks, galaxy centers, and ULIRGs are mostly 
captured by a surface density-dependent factor of the 
form /sB c>c S^jg^j, where Etotai refers to the combined 
gas plus stellar surface density on kpc scales. 

Present constraints remain scarce, but we make an ef- 
fort to present them in Fig. [121 The data corresponds 
to the kpc-scale dust-based measurements in nearby 
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disks by iSandstrom et al.l ()2012D . as well as the over- 
lap be tween the ULIRG samples bv lDownes fc Solomonl 
(|1998l fro m which we take d ynam ical masses and 
aco) and iPapadopoulos et aD (|2012L from which we 
take aco estimates). In this latter case, aco is de- 
rived from one-coin ponent (similar to the results by 
iDownes fc SolomorJ ) or two-component multi-transition 
fits (which include contributions from a dense phase). 
The dynamical surface density is do minated by the stel- 
lar component, even in ULIRGs (jPownes fc Solomonl 
IT99&) . Informed by the theoretical arguments lead- 
ing to Eg. [TBI and by the results of detailed mod- 
eling ()Shettv et al.l 120 lit) , we plot aco oc nor- 
malizing to our recommended Galactic aco value at 
Stotai = 100 M0 pc~^. Obviously this correction should 
not extend to surface densities below those of resolved 
self-gravitating GMCs. 

Given the large uncertainties and the small dynamic 
range of the aco measurements this simple prescrip- 
tion seems to reproduce the trends present in the data 
reasonably well, particularly for the results of one- 
component models for the ULIRGs (which we consider 
most mature). The observations may be fit with a 
smaller 7 although wi th considerable uncertainty (e.g., 
ISandstrom et al.ll2012[ ) , which leads us to prefer the the- 
oretically motivated 7 ~ 0.5. Density increases in the 
self-gravitating molecular material with respect to the 
Milky Way average GMC properties will drive the aco 
points up, while increases in temperature will drive them 
down. The sample spans a factor of ~ 2 in Tdust, which 
should be a reasonable proxy for gas temperature in the 
ULIRGs. Although we have searched for the signature 
of temperature effects in t he data, we see no di scernible 
correlation with Tdust (e.g.. lMagnelli et al.ll2012[ ). Likely, 
the sample lacks the necessary dynamic range to pull 
those effects out of the data. Possibly, as previously dis- 
cussed, the lack of a temperature correlation could be in 
part due to cancelations between the opposite effects the 
density of self-gravitating clouds and their temperature 
have on Xqq- Thus, as a tentative first step for a simple 
conversion factor prescription, we suggest using 



aco ~ 2.9 exp 



-0.4 



^ ^GMC 



■^total 



100 M© pc- 



(31) 



in M0 (K km s^^ pc^)^^, with 7 sa 0.5 for Stotai > 
100 M0 pc^^ and 7 = otherwise. Note that we still 
expect a fair dispersion around this average prescription, 
representing the variation in local parameters such as 
temperature or Sqmc- 

There has been an exciting range of theoretical and 
numerical developments on calculations of Xqq in the 
last few years. The coupling of high resolution hydro- 
dynamical simulations including chemistry and radiative 
transfer, with increasingly sophisticated theoretical mod- 
eling of photodissociation regions and molecular clouds, 
and galaxy scale simulations offers an exciting avenue 
of progress. Numerically derived c alibrations, such as 
those obt ained on sma l l scal es by iGlover fc Mac Lowl 
JlOllLor 'Shettv et al.l (|2011[ ) and on large scales by 
Narayanan et al. (2012), show much promise. Such sim- 
ulations are likely to become increasingly reliable as the 
modeling is able to better incorporate and couple the 
physics, kinematics, and radiative transfer on the small 
and large scales. Grounded on observations, simulations 
may offer the ultimate way to calibrate the CO-to-Il2 
conversion factor in a variety of environments. 
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